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ABSTRACT 


While  deuteron  stripping  theory  has  very  successfully 

explained  experimental  results  of  most  (d,n)  and  (d,p) 

» 

reactions,  there  are  some  instances  of  marked  disagreement  with 
the  simpler  forms  of  the  theory.  In  particular,  backward  peaks 
in  the  angular  distribution  of  the  reaction  products  sometimes 
appear;  an  explanation  of  these  based  on  an  exchange  effect 
called  heavy-particle  stripping  has  been  put  forward  by  Owen  amd 
JJadansky.  This  explanation  is  sound  in  principle,  but  involves 
an  over-abundance  of  adjustable  parameters  and  ancillary 
assumptions  which  make  detailed  verification  difficult  and 
ambiguous . 

The  present  work  contains  (l)  an  exposition  of  some  details 
of  the  basic  theory,  designed  to  make  what  follows  intelligible; 
(2)  a  critical  examination  of  work  published  by  the  proponents 
of  heavy-particle  stripping  concerning  the  reaction  most 
extensively  studied,  B;/  (d,n)C/2*4.43  Uev  {T)0IX,  including  the 
recent  extension  of  this  work  to  neutron-gamma  angular 
correlations;  and  (3)  an  account  of  an  attempted  alternative 
explanation  of  the  data  by  a  new  approximate  method  of 
int inducing  an  effect  of  Coulomb  forces  into  (d,n)  stripping 
theory;  this  attempt  does  not  succeed  in  explaining  the  Bli  (d,n) 
data,  but  the  method  is  of  3ome  interest  in  itself.  Other 
possible  explanations  are  also  mentioned. 
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Chapter  1  GENERAL  IlCT^ODUCTIOII 


A-  Description  of  the  Stripping  Process 

To  help  visualise  the  nature  of  reaction  mechanisms 
discussed  in  this  thesis,  a  brief  description  of  the  essential 
features  of  the  deuteron  stripping  proce33  is  presented  here  in 
qualitative  and  pictorial  terms.  Ho re  extensive  discussion 
along  this  line  can  be  found  in  the  first  chapter  of  Butler’s 
monograph  (l). 

The  usual  description  of  nuclear  reactions  in  terms  of  the 
compound  nucleus  hypothesis,  originally  due  to  ITiels  Bohr,  has 
very  successfully  described  most  but  not  all  nuclear  reactions. 
It  pictures  the  reaction  as  occurring  in  two  stages.  The  first 
stage  is  the  amalgamation  of  the  whole  incident  projectile  into 
the  target  nucleus  to  form  a  compound  nucleus,  not  directly 
observed,  and  generally  in  a  highly  excited  state.  After  a  time 
very  long  compared  with  the  time  it  would  have  taken  the 
projectile  3imply  to  traverse  target  dimensions  had  there  been 
no  interaction,  the  second  stage  of  the  reaction,  namely  the 
break-up  or  decay  of  the  compound  nucleus  into  outgoing  particle 
and  residual  nucleus,  takes  place  according  to  laws  independent 
of  the  first  stage.  The  compound  nucleus  retains  no  memory  of 
its  particular  mode  of  formation,  the  energy  and  momentum  of  the 
incoming  projectile  being  very  quickly  shared  with  all  nuclear 
constituents.  It3  decay  is  largely  akin  to  evaporation  of  a 
molecule  from  a  liquid  drop,  governed  entirely  by  statistical 
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laws.  There  is  thus  no  reason  why  reaction  products  should  come 
out  in  a  direction  close  to  the  incident  beam  direction,  and 
indeed  a  more  thorough  investigation  predicts  symmetry  between 
the  forward  and  backward  hemispheres  to  characterize  the  angular 
distribution  of  reaction  pruducts. 

Entirely  different  are  the  more  recently  studied  processes 
of  direct  nuclear  reaction,  of  which  deuteron  stripping  is  the 
best-known  instance.  There  is  no  compound  nucleus.  The  projectile 
usually  makes  a  grazing  collision  with  the  target  —  otherwise 
absorption  into  a  compound  nucleus  is  more  likely  to  occur  —  and 
interacts  only  during  the  short  time  it  takes  to  pass  the  target 
at  something  like  its  initial  speed.  As  a  result  the  outgoing 
particle  is  the  projectile  itself,  plus  or  minus  a  few  things 
lost  to  or  picked  up  from  the  target,  and  its  direction  i3  likely 
to  be  close  to  the  incident  beam  direction. 

The  deuteron  is  a  rather  loosely  bound  composite  structure 
in  which  the  constituent  neutron  and  proton  can  sometimes  be 

t 

found  separated  by  relatively  large  distances.  Because  of  this  a 
deuteron  incident  upon  a  nuclear  target  can  with  appreciable 
probability  make  such  a  collision  that  only  one  of  its  two 
constituent  nucleons  actually  hits  the  target  and  is  absorbed, 
while  the  other  misses  it  entirely  and  sails  on,  more  or  less 
undisturbed,  with  whatever  momentum  direction  it  had  at  the 
instant  of  collision,  to  be  counted  as  a  reaction  product, 

(Bigure  l)»  It  is  already  obvious  from  this  simple  description 
that  there  should  be  many  more  such  reaction  products  emitted  in 
directions  not  too  different  from  the  original  beam  direction 


. 
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Figure  1  Feuteron  Stripping 


Figure  2  The  Ilomentum  Balance  in  Feuteron  Stripping 


Figure  5  Heavy -Earticle  stripping 


than  in  the  "backward  direction.  This  will  be  confirmed  by 
detailed  evaluation  of  the  stripping  cross-section  in  the  next 
section. 

Whether  or  not  the  nucleon  which  hits  the  target  is 
absorbed  will  be  determined  to  a  large  extent  by  selection  rules 
corresponding  to  conservation  of  angular  momentum  and  of  parity. 
A  classical  interpretation  of  this  angular  momentum  immediately 
yields  more  information  about  the  angular  distribution.  The 
change  in  angular  momentum  between  the  initial  and  final  nuclei 
must  be  provided  by  the  total  angular  momentum  with  which  the 
Ccaptured  nucleon  comes  into  the  target;  this  in  turn  must  be 
compounded  of  that  nucleon's  spin,  of  magnitude  -J-,  and  of  its 
orbital  angular  momentum  with  respect  to  the  target  nucleus. 
This  fact  restricts  ^  to  very  fev/  possible  values;  conservation 
of  parity  restricts  it  further,  to  odd  values  only  if  the 
initial  'and  final  nuclei  have  different  parity,  to  even  values 
only  if  they  have  the  same  parity.  The  shape  of  the  angular 
distribution  depends  very  critically,  on  these  fev/  allowed  values 
of  yfc  —  in  practice  on  the  smallest  of  them  —  as  we  can  see 
classically:  7^  must  equal  kb,  where  k  is  the  relative  momentum 

of  captured  nucleon  approaching  target,  and  b  is  the  impact 
parameter,  v/hich  can  not  exceed  the  nuclear  radius  R  if  the 
nucleon  is  to  be  absorbed.  Thus 


Uow  the  momentum  k  is  provided  in  part  by  the  deuteron’s 
external  motion,  in  part  by  the  internal  motion  of  the  nucleon 
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while  "bound  in  the  deuteron,  so  that  the  probability  of  such 
momentum  being  available  decreases  sharply  with  k.  We  have  a 
fairly  good  idea  therefore  of  the  shape  of  the  cros3-section 
curve  as  a  function  of  k:  if  ^  —  0  is  allowed,  a  peak  occurs  at 
k  —  0  with  smaller  peaks  for  higher  allowed  fc  values;  if  fc  —  0 
is  not  allowed,  a  peak  occurs  at  k  ^  l/R  times  the  smallest  (c 
value  allowed,  with  further  smaller  peaks  for  higher  . 

But  it  is  easy  to  see  from  the  internal  momentum  balance  in 
the  deuteron’s  own  center  of  mass  that,  if  S  is  the  angle 
between  incident  beam  and  outgoing  reaction  product,  small  (9 
corresponds  to  small  k  and  large  S  to  large  k  (Figure  2):  for 
if,  as  is  most  probable,  the  outgoing  nucleon  stays  close  to  the 
deuteron  beam  direction  (small  d  ),  then  k  gets  little 
contribution  from  the  deuteron's  internal  motion;  whereas  if  the 
outgoing  nucleon  is  to  turn  around  completely  (large  O  ),  the 
captured  nucleon  must  get  a  big  push  forward  (hertce  large  k)  to 
conserve  momentum  in  the  deuteron's  center  of  mass.  Thus  the 
qualitative  description  of  the  cross-section  curve  as  a  function 
of  k  given  above,  also  describes  the  angular  distribution  of 
reaction  products  as  a  function  of  9  .  The  calculation  outlined 
below,  as  well  as  other  more  elaborate  ones,  confirm  these 
general  features. 


% 
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S-  The  Born  Anrroximation  Derivation 

Similar  formulas  for  the  angular  distribution  have  been 
obtained  by  Butler  (2),  using  an  elaborate  mathematical 
technique  to  match  wave  functions  and  their  derivatives  at  the 
nuclear  boundary,  and  by  Bhatia  et  al.  (5),  using  the  Born 
approximation.  Later  and  more  detailed  investigations  have 
established  that  Butler’s  approach,  while  sounder  in  principle, 
is  in  fact  equivalent  to  a  Born  approximation;  that  while  a 
Born  approximation  is  not  really  valid  here  in  the  usual  sense, 
there  are  reasons  which  condone  its  use;  and  that  either 
formula  accounts  strikingly  well  for  most  experimental  results, 
provided  one  is  allowed  to  treat  the  nuclear  radius  R  as  an 
adjustable  parameter  to  compensate  for  the  approximations  made. 
These  facts  will  be  taken  for  granted  without  discussion  in 
this  thesis. 

The  derivation  of  Bhatia  et  al.  is  much  simpler  than  any 
other,  leads  to  a  formula  which  can  be  conveniently  interpreted^ 
and  is  the  one  upon  which  the  work  reviewed  in  Chapter  2  has 
been  based.  The  outline  presented  here  follows  the  useful 

reviev/  article  by  Iluby  (4),  as  well  as  the  original  paper  of 

/ 

Bhatia  et  al.  (3). 

The  basic  approximations  common  to  all  treatments  must  be 
clearly  understood.  One  neglects  entirely  nuclear  interactions 
between  the  deuteron  and  the  target  nucleus  before  stripping, 
as  well  as  between  the  outgoing  nucleon  and  the  residual 
nucleus,  and  between  the  neutron  and  the  proton  from  the 


. 
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deuteron  after  stripping.  All  Coulomb  effects  are  likewise 
neglected  in  the  basic  treatment;  their  subsequent  introduction 
is  discussed  in  Chapter  3.  Under  these  assumptions,  products  of 

f 

plane  waves  and  appropriate  internal  wave  functions  can  be  used 

for  the  initial  and  final  wave  functions  to  be  inserted  into 

the  born  formula.  The  basic  idea  of  this  treatment  —  and  of 

► 

others  —  is  that  the  deuteron  is  broken  up  instantaneously  and 
its  momentum  balance  at  that  instant  determines  subsequent 
happenings:  the  corresponding  momentum  of  the  outgoing  nucleon 
is  the  momentum  actually  observed,  which  defines  9  ;  and  in 
effect,  the  captured  nucleon  is  then  pictured  as  approaching 
the  target  nucleus  ns  a  plane  wave  with  the  momentum  it  obtains 
from  that  break-up. 

The  formal  derivation  is  presented  in  this  section,  and 
the  significance  of  the  various  factors  is  discussed  in  the 
next.  Consider  the  reaction  A(d,b)B  (Figure  l),  whereby 
deuteron  bombardment  of  an  initial  nucleus  A,  with  spin 
leaves  the  residual  nucleus  B  in  a  particular  state  of  spin  J^. 
The  subscript  c  will  indicate  quantities  belonging  to  the 
nucleon  c  v.Thich  is  captured  (e.g.  the  proton  in  a  (d,n) 
reaction) • 

According  to  the  Born  approximation,  the  cros3-3ection  for 
the  reaction  to  proceed  from  initial  magnetic  substates  of  A 
andyU^  of  the  deuteron  to  final  substates  lEg  of  B  and  ^  of  the 
outgoing  nucleon  is  found  to  be 


B1 
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<• 


where 


■]%<£&*)  gW* 


Thi3  expression  gives  tne  cros3-secf ion  in  the  center- of -mass 
coordinates  provided  the  mass  symbols,  are  given  the  following 
meanings : 


*  _  . 
md^mA 


n^mg 

■■  '  ■  ■■  • 


^A 

mc+IaA 


is  the  initial  wave  number  of  relative  motion  of  deuteron 
with  respect  to  A,  so  that  the  total  kinetic  energy  before 
collision  in  the  center- of -mass  system  is 

Ti2  kf 
2  md*  ‘ 

Similarly,  k^  is  the  wave  number  of  relative  motion  of  b  with 
respect  to  B  after  collision,  the  corresponding  kinetic  energy 
being  given  by 

h2  kf 

in  the  center-of-mass  frame.  X*.  X. 

}  A 

functions  of  the  deuteron,  target  nucleus  A  (ground  state),  and 
residual  nucleus  B,  respectively.  V(rc,cxr,  £)  is  the  interaction 
potential  between  the  nucleons  in  A,  whose  coordinates  are 
collectively  denoted  by  t  ,  and  the  captured  nucleon  c.  And  rc, 

C Tc,  r^j,  (T-^j  etc.  are  respectively  the  position  and  spin 

coordinates  of  the  particle  designated  by  the  subscript. 


Z 


tii 


B 


are  internal  wave 
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T  • 


V/e  pick  out  first  from  B2  the  integral  over  d  £  , 

/ 

~X  H  ^ 9  rc »  }  )  v  >  X)  (  /  )  d  /  » 


which  is  a  function  of  rc  and  of  <7~c,  and  assume  that  the 
capture  takes  place  at  the  radius  rc  —  B,  so  that  we  write 

cT (rc  -  R) 


T*  — 


~ n -  T" 

R^ 


where  T"  is  a  function  of  Oc  *  c ,  and  (T z  hut  is  independent 
of  rc,  and  may  thus  he  expanded  as  follows: 


iji  ii  = 


The 


Dirac  'brackets  <.  11/  are  numerical  coefficients 


<%  I  Vl  %»-£>“c-Ac>  Yl  (*c  fc)  JA^) 


giving  the  projection  of  T”  on  the  product  function  Xj  c  • 

The  Y^c  are  the  usual  spherical  harmonics  and  the  'Xmc  are 
eigenfunctions  of  a  spin  -J-  particle  with  components  //c  along 
the  axis  of  quantization. 

Substituting  hack  one  finds 


T  —  *\UB  I  V  I  J A’  >mc  y^c’  b  /  dy> 


B3 


^ c  ^c/^c 
where 


bid y s  J 

*  X  ^  drh  a'ir ^  dcoc 


B4 


The  symbol  R  denotes  the  vector  whose  spherical  polar 
coordinates  are  rc  —  R,  G  c,  <^c.  The  element  of  solid  angle 
sin  &  c  d^c  d  f  q  is  written  d^c. 
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The  expression  discussed  above  for  | T / 2  should  really  be 


summed  over  final  and  averaged  over  initial  magnetic  quantum 
numbers  ^  dJ  *-A>  an^  -4b*  3ince  these  quantities  are 

experimentally  unspecified:  the  various  final  substates  all 
contribute  incoherently,  and  the  incident  beam  and  target 
nuclei  are  unpolarized.  Hence 


/yfr/K 


The  3  and  (2Ja  +  l)  are  the  number  of  magnetic  substates  of  the 
deuteron  and  target  A  respectively. 

neglecting  the  d-state  admixture  in  the  deuteron  v/ave 
function,  one  use3  for  the  product  of  a  triplet  spin  factor 
and  an  s-state  (spherically  symmetric)  orbital  factor: 


Equation  B4 


where 


and 
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is  the  Clebsch-Gordan  coefficient  corresponding  to  the 
composition  x  i  of  the  proton  and  neutron  angular  momenta 
with  componentsy6/c,  into  the  deuteron  triplet  state  v/ith 
component Hence 


|rr 


Ma  /  [l  V‘)  {/'Alfa)  fa) 


The  summations  over  >  and  /ut  are  considerably 

simplified  by  properties  of  the  Clebsch-Gordan  coefficients: 
^ujJhj  Ud  /  ”  0  unless  /ty*  Therefore  for  a  given 

term  in  the  outside  summation,  there  is  only  one  non-zero 

term  in  the  coherent  summation,  which  may  thus  be  performed 

incoherently.  Moreover  one  really  sums  over  only  two  of  the 

/ 
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three  quantum  numbers  and since  they  are  so 

related.  Expressing  it  as  a  3um  over  and  leads  to 


IrlA 


i 


»,A)  fafa 


are  zero 


/  \  .  (fa  A*  i  fad) 

ept  /  [  fa  ~  \  I  -fa  =  1.  and  <fa,  ■  o)  — 

fa-i,i  l  fa*) 


How  all  Clebsch-Gordan  coefficients 
exc 


fap  .  Inserting  these  values  gives 


Hollowing  13hatia  et  al.,  the  summations  over  uc»  11^,  1%, 
can  now  be  carried  out  with  the  aid  of  the  group-theoretical 
properties  of  the  wave  functions  concerned  in  the  definition  of 
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the  Dirac  bracket  expressions 
result  is 


<^tB  V  i  ^c  ,mc  >/^c^.  ’  • 


The 


where 


/l 


% 


'sM8\v\Hu</*tyi 


V<c  /*<. 

t  u  “b 


is  usually  kept  as  an  adjustable  parameter,  since  its  honest 
evaluation  would  require  more  detailed  knowledge  of  nuclear 
structure  than  is  likely  to  be  available  for  some  time. 

Tfi 

Finally  g-fc  can  be  evaluated  from  its  definition.  Equation 

A(, 

B5.  Yftien  performing  the  £b  integration,  rc  r?  E  -  (E,  <9C,  /  c )  is 
constant,  so  that  the  integration  over  r^  can  be  replaced  by  an 
integration  over  the  relative  coordinate  ft  -  r^  -  E  : 

=  5 yf(6. V % (Ve dr ** 

-  o-(Kt)  i 'yphtje-  -  </*'. 


where  Kj,  -  -  -g-k^  represents  that  part  of  the  outgoing 

nucleon fs  momentum  which  is  obtained  from  its  internal  motion 
in  the  deuteron;  whereas  k  —  k^  -  k^  (or  more  accurately 

k  =  k^  -  kb  finite  mas3  effects  are  considered)  is  the 

“  —  — 

momentum  with  which  the  captured  nucleon  c  approaches  and 


15 


enters  the  target  nucleus  A  (figure  2).  The  symbol  G (K^ ) 
denotes 


0  w  *  *  ,,  *k'l 

Finally,  J  Jjp  c  (v^Tc)£  <J(*{  is  readily  evaluated  if 

the  axis  of  quantization  for  determining  is  taken  in  the 
direction  of  k.  Starting  with  the  well-known  expansion  of  a 
plane  wave  into  spherical  harmonics. 


where 


is  the  spherical  Bessel  function,  one  obtains 


y  >>77 -(ll-n) JjJkR)  cC£/ 


Hence , 


G  V Q  hr  -v)  j  a  (kn) 


and  substituting  back  into  Equations  B6  and  Bl, 


<r 


>  /-w*  _  y»y>.  k_ 

+vr  *  z  +  1  1  ~  4, 


. 
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G-  Discussion  of  the  Cross-Section  Formula 


Equation  B9,  the  final  formula  for  the  stripping  cross 
section,  is  conveniently  expressed  as 


<r  =  TT  m  ^  Lftt)  Pf< 


Cl 


v/liere 


W*i)  =  M~3Gl(f6) 

,'k-R 


C2 


L/c  2  / 1 \  ‘AtU  e/cvcj  =  4^/2^+/; jjJkt)  C3 


and 


6 


TT  <  »/.  h  Ah 

t¥  ~£  (^)diH) 
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These  three  factors  have  patent  physical  significance. 

The  first  factor  TTfe fc)  describes  the  break-up  of  the 
deuteron.  It  is  the  Bourier  transform  of  the  deuteron's 
internal  wave  function  and  thus  gives  the  probability  for  the 
outgoing  nucleon  b  to  obtain  from  the  deuteron’s  internal 
motion  the  momentum  K-^  which,  added  to  half  of  the  deuteron'3 
external  momentum  k^,  produces  the  observed  outgoing  momentum 
k^. (Bigure  2).  JJ (K^)  decreases  smoothly  and  rapidly  with 
increasing  0  ,  and  is  not  the  factor  mainly  responsible  for  the 
typical  peaked  shape  of  the  cross-3ection  curve  described 
qualitatively  in  Section  A. 

L j;  (k)  is  the  centrifugal  barrier  factor  for  a  nucleon  of 
momentum  k  to  reach  the  target  boundary  with  angular  momentum 
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as  discussed  below.  And  P^  is  the  probability  that  if  such 
a  nucleon  does  reach  the  target  boundary  it  will  be  absorbed* 

P p  is  independent  of  k  and  of  9  . 

The  characteristic  peaks  of  the  angular  dependence  are 
determined  mainly  by  the  centrifugal  barrier  factor  Iy  (k).  The 
significance  of  this  factor  is  clear  from  the  nature  of  the 

expansion  of  a  plane  wave  into  spherical  harmonics  used  above 

n  -y,  'L  ft 

(B8)  to  evaluate  the  integral  j  ^  C(%(fc)e  .  Prom 

general  principles  of  quantum  theory,  the  modulus  squared  of 

/  / - 7 

the  coefficient  i  c  j/4^(2/c  -f  1)  (kR)  of  y/  (  &c  <t>c  )  in  this 
expansion  then  gives  the  probability  that  the  system  described 
by  the  plane  wave  is  in  the  state  Y y  (  9C  <pc  )  at  R.  But  this 
modulus  squared  is  precisely 

p _ __ _ _ ^  1 2^ 

i  ‘^(24  +17  (Id i)  =  4  77  (2^.  +  1)  (kK)-  Lyt(k)  . 

Thus  the  nucleon  c  may  be  pictured  as  issuing  from  the  deuteron 
as  a  plane  wave,  which  is  equivalent  to  this  sum  of  spherical 
waves  of  which  only  a  few  are  effectively  absorbed,  since  the 
selection  rules  make  P^  —  0  for  all  but  a  few  values  of  vf . 

The  properties  of  the  nucleus  exert  their  influence  upon 
the  probability  of  the  reaction  taking  place,  through  these 
selection  rules  on  the  Tp  .  The  nucleus  having  thus  declared 
which  of  the  states  "^Tc  it  i3  willing  to  accept  and  to  what 
extent,  the  factor  (k)  determines  how  much  of  each  such 
3tate  is  available  in  e^*-,  the  plane  wave  representing  the 
approach  of  the  nucleon  c.  The  factor  ff  (K^ )  represents  the 
influence  of  the  dynamics  of  the  deuteron  break-up  in 


. 
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determining  the  intensity  of  the  plane  wave  e-  -  as  a  function 
of  the  relative  momentum  k. 

In  a  way  it  is  rather  surprising  that  such  an  interpretation 
should  he  a  valid  description,  for  the  picture  of  a  plane  wave 
of  nucleons  c  is  even  more  virtual  and  subjective  (in  the  sense 
of  being  a  device  for  understanding,  created  by  the  physicist’s 
mind)  than  most  of  the  other  constructs  of  physics.  In  the 
first  plaee  these  three  events  —  break-up  of  the  deuteron, 
approach  of  c  towards  A  in  the  form  of  a  plane  wave  with 
momentum  k,  and  absorption  of  c  by  A  —  are  separated  only  for 
convenience  of  description:  they  are  not  sequential  in  time. 

This  makes  the  intermediate  plane  wave  considerably  less  real 
than  the  compound  nucleus  in  a  Bohr-type  nuclear  reaction  for 
instance,  which,  although  not  observed,  must  be  postulated  to 
exist  for  a  specified  length  of  time.  Moreover  the  nucleon  c 
must  enter  the  nucleus  with  negative  energy  for  a  bound  state 
to  result,  so  that  it  is  only  by  analogy  that  one  can  think  of 
this  plane  wave  as  a  solution  of  the  free-particle  Schrddinger 
equation,  since  the  momentum  k  does  not  correspond  to  an 

o  p 

energy  in  the  usual  way. 

2mc 

On  the  other  hand  it  should  be  emphasized  that  the 
experimental  evidence  confirms  this  formula  very  clearly.  The 
formula  predicts  a  very  special  and  recognizable  angular 
dependence,  governed  by  these  hj  (k)  coefficients  in  the 
plane  wave  expansion,  and  the  fact  that  the  formula,  admittedly 
approximate  and  far  from  unique,  nevertheless  fits  the 
experimental  data  very  well,  is  an  encouragement  to  pursue  this 
analogy  in  attempting  to  extend  the  theory.  It  should  be  borne 


■ 
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in  mind  that  the  angular  distributions  are  determined  primarily 
by  the  momentum  balance,  and  the  momentum  balance  is  preserved 
artfully  in  this  perturbation  approach.  Ueutron-gamma  angular 
correlations  discussed  in  the  next  section  also  support  this 
picture. 

It  is  interesting  that  the  Butler  formula  obtained  by 
smooth  fitting  of  wave  functions  at  the  nuclear  boundary,  can 
also  be  cast  into  the  form  Cl, 

<r  =  racj  LJt  an  p£e 

Yc 

and  30  shares  in  the  above  interpretation  to  a  certain  extent, 

/ 

if  less  directly.  The  factor  TJ  (K^)  is  unchanged.  Instead  of 
involving  simply  a  spherical  Bessel  function  j  j  (kR) ,  the 
factor  L  £  (k)  in  the  Butler  formula  involves  the  Wronskian  of 

t 

j  p  (kR)  ylth  a  function  independent  of  angle,  so  that  it 
contains  also  a  correction  term  proportional  to  the  derivative 
of  j  £  (kR).  This  correction  term,  which  is  in  fact  dominated  by 
the  first  term  in  all  cases  where  the  basic  stripping 
approximations  are  justified,  can  be  neglected  in  most 
practical  applications,  and  its  effect  simulated  by  using  a 
radius  artificially  larger  by  about  one  fermi  in  the  Born 
approximation  formula  than  in  the  Butler  formula.  The  effective 
radius  in  both  theories  appears  more  a3  an  adjustable  parameter 
than  as  a  physically  meaningful  quantity. 

The  Butler  derivation  gives  for  ,  in  the  case  of  a 
bound  final  state, 

n-R2m^kb  (2  Jb  t  1) 

\  '  —2  ^  (2J4  +  l)2(2/e  +  1)  \ 


C5 
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where  is  the  probability  density  that  the  state  of  B  can 

yield  c  at  the  surface  with  the  orbital  angular  momentum  /  , 
leaving  A  in  its  ground  state,  llore  precisely,  let  the 
normalized  wave  function  js[1  (rc,  <TC,  'f  )  in  the  region  rc  ^  H 
(where  B  breaks  up  into  the  non-interacting  constituents  A  and 
c)  be  expanded  as 


'^V//  ^ 

where  the  functions  Uy.  *.  (  )  describe  the  states  of  A  (with 

angular  momentum  JA ,  liA),  and  .  .  (rr)  is  a  free- 

particle  radial  wave  function  for  c.  Then  P  is  given  by 


!  „  f  f 
1  n  *  ,  A 

*  4  -V  »-  '  c 


r*; 


Ha,  *S/4 

where  t  =  0  designates  the  ground  state  of  A. 

On  the  other  hand,  the  Butler  result  for  the  case  of  a 
virtual  (instead  of  bound)  level  of  B  is  given  by 


TT  kh  (xTgirl) 


K. 


f  kd  C^OzCz/c*/) 

Here  ^  is  the  reduced  width  of  the  level,  related  to  its 

partial  width  fj  (for  decay  into  c  and  the  ground  state  of  A, 

C  * 

with  angular  momentum  ^  )  by 


C6 


y/  =  i{xU2<wr  I}. 

One  can,  by  analogy,  formally  define  a  "reduced  width"  for  the 
bound  levels  also, 


o 

I  /C. 


„„ 
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.  (i) 

which  makes  Equation  C 6  coincide  with  C5.  (  ly  is  a  spherical 
Harikel  function  of  the  first  kind,  and  k^  is  such  that  ^c^c  i3 
the  energy  of  the  virtual  level;  kc  is  not  the  same  as  k,  and 
does  not  depend  on  0.)  The  Born  approximation  formula  C4  for 
P y  can  he  reconciled  with  the  Butler  form  C5  by  formal 
manipulations* 

To  close  this  discussion  let  us  note  from  Equation  B9  that 
the  contributions  of  the  various  values  of  Ji  allowed  by 
conservation  of  angular  momentum,  i*e. 


C7 


are  to  be  summed  incoherently,  and  that  in  practice  all  but  one 
of  these  values  can  be  neglected;  the  one  important  value  is 
usually  the  smallest  allowed,  except  in  some  interesting  cases 
where  this  is  excluded  by  nuclear  properties.  These  cases  are 
usually  explicable  in  terms  of  the  shell  model. 


p-  Angular  Correlations 

The  theory  has  been  extended,  by  Satchler  (5)  and  by 
others^  to  describe  the  angular  correlations  between  the 
direction  of  emission  of  nucleon  b  and  that  of  the  subsequent 
gamma  ray  of  de-excitation,  in  cases  where  B  is  left  in  an 
excited  state.  An  angular  correlation  function  gives  3imply  a 
coincidence  counting  rate  as  a  function  of  the  two  directions: 
hoY/  many  gamma  rays  are  recorded  in  one  direction  in 
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coincidence  with  particles  b  in  another  direction. 

It  i3  an  interesting  confirmation  of  the  picture  presented 
in  the  previous  section  that  the  angular  correlations  can  be 
accounted  for  on  the  following  basis.  Particle  c  is  pictured  as 
entering  nucleus  A  as  a  plane  Y/ave  of  momentum  k,  which 
momentum  is  related  to  the  direction  of  emission  of  b;  the 
subsequent  angular  distribution  of  gamma  ray3  with  respect  to 
the  direction  k  is  calculated  exactly  on  the  same  basis  as  if 
one  were  studying  the  reaction  c  «f  A  — *  (B  )  — *■ >  B  ■+•  Y  by  means 
of  bombardment  of  A  by  a  beam  of  particles  c  incident  in  the 
direction  k. 

It  follows  that  the  angular  correlation  curves  should  be 
symmetrical  with  respect  to  this  direction  k,  which  is  the 
recoil  direction.  This  general  property  of  stripping  angular 
correlations  is  usually  confirmed  by  experiment  v/henever  the 
angular  distributions  show  characteristic  stripping  patterns. 

*  The  distribution  of  gamma  ray  directions  following  such 

i* 

bombardment  is  given  quite  simply,  in  cases  of  pure 
multipolarity  and  short  enough  lifetime  (  T  ^  10  U  sec.)  to 
avoid  complications,  in  terms  of  functions  v/hich  will  now  be 
defined.  This  outline  follows  mainly  an  unpublished  report  by 
Edwards  (6)  which  gives  in  les3  up-to-date  but  somewhat  more 
explicit  form  an  account  of  his  published  work  (7). 

The  pure  multipole  fields  are  a  complete  set  of  solutions 
of  the  free-space  Eaxwell  equations,  characterized  by  their 
parity  properties.  Instead  of  going  through  the  general 
solution  of  this  problem  involving  the  vector  spherical 
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harmonics,  vie  need  to  consider  only  a  more  restricted  and 
easily  obtained  set  of  functions  capable  of  describing  the 
simpler  fields  prevailing  in  the  so-called  "far  field"  or  “wave 
zone"  region,  i.e.  far  from  the  nucleus  which  is  the  source  of 
the  field.  In  that  region,  the  fields  £  and  are  almost 
transverse  and  the  propagation  vector,  or  the  vector  G  giving 
the  linear  momentum  density  of  the  field, 


£  x  H 

‘irrc. 

which  is  proportional  to  it,  i3  almost  parallel  to  the  radius 
vector  r. 

We  can  not,  however,  make  the  simplifying  assumption  that 
^  and  ft  are  both  perpendicular  to  r,  because  we  are 
specifically  interested  in  describing  the  removal  of  a  certain 
quantity  of  angular  momentum  through  emission  of  a  photon,  and 
the  angular  momentum  density  of  the  field 


r  x  G 

477C 


of  course  vanishes  under  these  oversimplified  conditions. 

On  the  other  hand  we  can  assume  that  only  one  of  the  two 
vectors  £  and  ^  has  a  non- vanishing  component  in  the  r 
direction,  while  the  other  is  always  perpendicular  to  r.  Then 
since  &  and  fb  are  perpendicular  to  each  other  and  equal  in 
magnitude  in  the  far  field  region,  the  angular  dependence  is 
given  by  the  magnitude 


^  —  fi  * 

f7 r  i77  ' 
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of  the  Poynting  vector 


?  -  £  f 

We  can  substitute  into  D1  an  expression  for  either  or 
whichever  is  the  simpler. 

Let  us  construct  then  a  complete  set  of  vector  point 
functions  of  &  and  f  ,  everywhere  perpendicular  to  r.  Taking 
the  gradient 
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,  a  ,  p 
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r  S'n  0 


of  the  scalar  (ordinary)  spnerical  harmonics  h[9f)  is  the 
simplest  way  to  obtain  a  set  of  vector  point  functions  of  9 
and  ;  cross-multiplying  by  r  eliminates  the  r-dependence 
introduced  by  the  gradient  operation,  and  makes  the  functions 

f 

always  perpendicular  to  r: 


A  x  t7  — 


$ 


P 


-  07T 
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d&  v  u'n  9 


For  convenience  the  set  of  functions  XLn{fif)  is  defined  as 

Xu,  (9J1-  2-  * 1  i  ffj  >  ' 


\iLO--ri )  ZLO.-H)'  C  - '/ 

r-r 

Because  the  operator  r  x  _  commutes  with  the  angular  momentum 
operator  of  which  the  YL/1  are  a  complete  orthogonal  3et  of 
eigenfunctions,  these  are  a  complete  orthogonal  set  of 

functions  perpendicular  to  r. 

To  discuss  the  case  where  the  vector  rfc  is  perpendicular 
to  r  in  the  far  field  region,  we  define 


' 
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It  then  follows  from  the  free-space  Maxwell  equation  that  the 
corresponding  £  is 


i 


£ ) 

•  /  * 


How  the  coefficients  g  (r)  can  he  suitably  chosen  so  that 
the  H  satisfies  the  free-space  Harwell  equation.  koreover, 
since  the  YA.  (£,/■)  have  parity  (-1)",  it  follows  that  ;*  (r) 

(  C 

as  defined  has  parity  (-1)'  and  (r)  has  parity  -(-l)~.  The 

— — - 

statements  of  this  paragraph  amount  to  a  definition  of  what  is 
known  as  an  electric  multipole  field  of  multipolarity  (L,Ii). 

On  the  other  hand  we  could  have  started  by  assuming  that 
the  vector  ~  instead  of  £*.  was  perpendicular  to  r  in  the  far 
field  region,  and  constructed  a  vector  point  function 

C  (I)  =  &  *> 


with  the  g  ,,(r)  suitably  chosen  to  satisfy  the  Harwell 
U  ^ 

equations  for  £  .  This  £c  A  (r)  has  parity  (-1)^  and  is 
accompanied  by  an  $£/r(r)  having  parity  -(-l)^.  This  set  of 

l)  and  £)  then  corresponds  to  the  definition  of  what 

is  known  as  the  pure  magnetic  multipole  field  of  multipolarity 


Having  obtained  such  expressions  for  the  electric  and 
magnetic  multipole  fields,  it  is  a  simple  matter  to  obtain  the 
angular  distribution  of  the  corresponding  gamma  rays,  through 
that  of  the  Poynting  vector  giving  the  energy  flow,  yor  the 
electric  multipole  field,  one  finds 

iXh 
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and  similarly  for  the  magnetic  multipole  fields: 

ur.|=4C‘=  writ w-x Lie.*). 

Thus  in  both  causes  the  angular  dependence  is  contained  in  the 
functions 

f[,„(  M 

=  uto  {  jlte  ^17 

which  are  called  the  classical  multipole  distribution  functions. 
Tables  of  these  are  available  (Arnold,  8). 

In  applying  this  theory  to  the  stripping  reaction  leading 
to  the  first  excited  state  of  B,  with  spin  Jg  ,  which 
subsequently  decays  to  the  ground  state  with  spin  (subscripts 
e  for  excited,  f  for  final  state)  by  a  one-step  pure  multipole 
transition  of  multipolarity  (Lfli),  it  will  be  necessary, 
however,  to  perform  coherent  summations  over  the  possible 

values  of  any  unobserved  intermediate  quantum  numbers,  namely 

* 

the  various  possible  combinations  of  magnetic  substates.  Hence 
one  must  U3e  the  appropriate  amplitude  expressions 
instead  of  the  probability  expressions  ^^(<9)  in  order  to 
obtain  the  coherent  effects  on  the  angular  distributions. 

In  particular  if  the  system  is  in  a  magnetic  substate  lle 
of  Je  9  the  amplitude  function  corresponding  to  its  decay  to  a 
particular  substate  of  i3  given  by  XLfJ  ( dj  <f>) ,  with 
II  a  -  1> ,  and  this  is  to  be  summed  over  the  various  possible 

1!£  with  weighting  factors  giving  the  relative  probability 

T  One  can  show  (46)  that  the  multipole  fields  as  defined  here 
do  carry  away  angular  momentum  (Xj,h)  • 
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amplitude  of  the  transitions  ii^  to  Up. 

Let  re(H),  [f  (R)  represent  the  radial  functions  of  the 
internal  coordinates  of  B  corresponding  to  excited  and  final 
states  respectively,  ,Mp)  represent  the 

corresponding  tota.1  angular  momentum  eigenfunctions,  and 
^(Ljh)  the  wave  function  for  the  photon  emitted.  The  internal 
wavefunction  of  the  excited  state  can  then  be  written 
/g  (R)  ^(Je,He)  and  that  of  the  final  state  (R)  ^,(  ,M^)  ^(L,ll) . 

The  angular  part  of  this  final  3tate  wave  function  is 
conveniently  expanded  in  terms  of  excited  state  angular 
functions : 


l  Gw  =2L  c%%L"  &  u'^'} 

£  ‘‘V 

r-^  *r  <- 

where  the  Clebsch-Gordan  coefficient  corresponding 

to  Je’  =  L  with  Mi  •=-  ly  -f  M.  The  relative  probability 

amplitudes  of  the  transitions  are  then 


''  nt  m 


and  the  amplitude  for  the  angular  distribution  resulting  from 
the  transition  Me  -vllp  becomes  proportional  to 

A  (Mo  i 9 ‘p)  ~  V  (04>) 

r\y~  I  '£/  t  /  —  * 


This  amplitude  must  be  multiplied  by  a  further  weighting 
factor,  namely  the  probability  amplitude  that  the  stripping 
process  will  leave  B  ^ in  substate  and  then  summed 
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coherently  over  the  possible  values  of  1^.  This  stripping 
amplitude  can  be  evaluated  for  instance  from  Equation  B3,  with 
1£ ^  =  U  ,  except  that  one  normally  neglects  the  sum  over  fc  , 
using  only  the  one  value  of  <c  which  is  known  from  the  angular 
distribution  to  be  predominant;  measurements  are  usually  taken 
with  the  b-counter  at  the  major  peak  of  the  angular 
distribution  corresponding  to  that  Ac  ,  so  that  this  procedure 
should  be  justified.  Final  expressions  must  as  usual  be  summed 
over  possible  values  of  the  final  and  averaged  over  those  of 
the  initial  unobserved  quantum  numbers,  so  that  the  angular 
correlation  function  is  proportional  to 

%  %,%) 

'  ‘H  *  *  f  U 

where  T  (0^)  denotes  the  stripping  amplitude  or  matrix  element, 

as  in  Equation  B2  with  lit  =■  K  .  (  refer  to  the  gamma-ray 

direction  and  9^  ,  written  9  in  other  sections,  to  that  of  b.  ) 


/ 

C 1 . 
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E-  The  Reaction 


x> 


11  fd.n)Ca^4.43  Uev 


The  basic  stripping  theory  presented  in  the  preceding 
sections  has  proved  extremely  successful  and  useful  in 
accounting  for  experimental  angular  distributions  of  many  (d,n) 
and  (d,p)  reactions,  not  only  under  conditions  where  its  basic 
assumptions  are  justified,  but  even  in  many  cases  where  they 
are  not.  There  are  several  reactions  however  for  which  the 
experimental  data  do  not  exhibit  the  usual  charact eristic 
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stripping  pattern.  Because  of  the  success  obtained  in  other 
cases,  attempts  have  been  made  to  account  for  the  anomalous 
results  by  modifications  and  extensions  within  the  general 
framework  of  deuteron  stripping  theory. 

One  of  these  reactions,  B 11  (d,n)C'""  leading  to  the  first 

i  x 

excited  state  of  C  at  4,43  Uev,  forms  a  major  subject  of  this 
thesis.  The  experimental  data  are  in  the  range  of  bombarding 
energies  from  ^  to  2  Hev* 

The  B/;  ground  state  has  spin/parity  3/2",  as  expected  on 
the  basis  of  the  shell  model.  The  ground  state  of  C  ,x  has  0  ", 
and  the  4.43  Kev  first  excited  state  2%  as  is  normal  for  even- 
even  nuclei.  (The  gamma  ray  of  de-excitation  is  pure  quadrupole 
as  it  should  be.)  Conservation  of  total  angular  momentum  in  the 
formation  of  the  2^3tate  thus  restricts  to  0  ic  4 

according  to  Equation  C7 ;  and  conservation  of  parity  further 
restricts  sfc  to  the  odd  values  ~  1  and  3. 

Prom  pure  stripping  theory,  the  angular  distribution 
expected  (Cl)  is  therefore  proportional  to  TF{K^)  j  a (kR) 
plus  some  minor  contribution  of  TT(lZ^)  j^(kH).  Por  \  Eev,  3ay, 
the  distribution  should  start  from  a  minimum  (seldom 

actually  seen)  at  <9  —  0,  rise  to  a  peak  around  10° to  15*, 
decrease  sharply  to  a  zero  around  105°,  and  then  rise  slightly 
towards  a  flat  secondary  maximum  around  175^,  thi3  secondary 
maximum  being  only  about  6/  as  high  as  the  forward  maximum. 

The  experimental  curves  actually  obtained,  e.g.  by  lleilson 
et  al  (9),  are  very  different.  They  show  a  maximum  at  9-0, 
decrease  only  slightly  to  a  minimum  whose  position  varies  from 


■ 
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about  6b°  at  Hev  to  about  90°  at  2  Hev;  then  ri3e  again  for 
larger  angles.  The  ordinate  at  the  minimum  i3  from  70%  to  50^ 
of  the  height  of  the  forward  peak,  instead  of  going  down  to 

"zero",  and  the  backward  maximum  is  from  80%  to  60 %  as  high  as 

o 

the  forward  at  the  highest  angles  reached,  about  145  . 

The  experimental  angular  distribution  for  a  bombarding 
energy  of  -J-  Hev  is  presented  in  Figure  4,  along  with  a 
stripping  curve  for  these  conditions.  The  data  at  the  other 
bombarding  energies  show  similar  disagreement  with  theory.  It 
wall  be  more  convenient  to  present  them  for  detailed  study  in 
a  different  form  (plotted  as  a  function  of  k  instead  of  0  )  in 
Chapter  3. 

Heutron-gamma  angular  correlations  for  this  reaction  have 
also  been  measured  by  Keilson  et  al  (9).  They  do  not  show 
syasnetry  about  the  recoil  direction  k,  as  they  should  on  the 
basis  of  pure  stripping  theory  (see  Section  D). 

A  reasonable  radius  for  the  nucleus  B 11  is  of  the  order  of 
3.5  fermis,  which  implies  that  the  deuteron,  or  at  least  its 
proton  constituent,  must  approach  to  a  radius  of  about  4.5 
fermi3  if  it  is  to  be  stripped  off  by  the  attractive  nuclear 
forces.  The  height  of  the  Coulomb  barrier  as  seen  by  the  proton 
at  this  ra-dius  would  be  of  the  order  of  1.6  Hev  even  for  fc^Q. 
For  Jl^-  1  the  centrifugal  barrier  increases  the  barrier  height 
to  4.1  Hev,  at  a  distance  of  4.5  fermis.  Hence  it  would  not  be 
unreasonable  to  expect  the  stripping  pattern  to  be  very 
considerably  modified  by  Coulomb  effects  at  the  energies 
considered.  On  the  other  hand,  there  are  other  reactions  where 
a  recognizable  stripping  pattern  persists  even  under  such 
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adverse  conditions* 

Recent  work  by  Wilkinson  (10)  lias  suggested  reasons  why- 
stripping  reactions  with  low  ^-values,  sey  from  -2  to  +  2  Hev, 
might  he  expected  to  exhibit  persistent  stripping  patterns  even 
at  low  energies.  The  reaction  B /  (d,n)C  x  has  a  high  Q-value 
of  9.50  for  this  level  and  so  is  not  in  that  category;  on  the 
contrary  the  Wilkinson  argument  would  predict  a  much  le33 
characteristic  stripping  pattern  for  such  a  high  Q- value. 

The  explanation  of  these  anomalous  distributions  by  heavy- 
particle  stripping  will  be  critically  reviewed  in  Chapter  2, 
and  the  possibility  of  other  explanations  investigated  in 
Chapter  3. 


r* 
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(ftia-oter  2  HEAVY -PART  I C LB  STRIPIIUG 

g-  Peacrlption  of  Heavy-rartlcle  Stripping 

The  unexpected  angular  distributions  obtained 
experimentally  for  B;/  (d,n)C /z  *(4.43  Hev  level)  have  been 
described  in  Section  E  above.  These  and  similar  results 
obtained  in  a  few  other  reactions  have  led  a  group  at  Johns 
Hopkins  University  to  postulate  an  effect  called  heavy-particle 
stripping.  They  have  published  several  papers  in  connection 
with  this  proposal  (ll  to  19,  7).  This  hypothesis  will  be 
described  qualitatively  in  this  section,  and  examined  more 
critically  in  the  remainder  of  this  chapter. 

The  essential  idea  of  the  ordinary  deuteron  stripping 
mechanism  in  a  (d,n)  reaction  is  that  the  outgoing  neutron 
comes  from  the  deuteron,  not  the  target  (Figure  1).  This 
neutron  bears  no  label  of  course,  but  it3  identity  is  established 
by  the  success  of  deuteron  stripping  theory  in  predicting  such 
data  as  angular  distributions.  But  suppose  now  the  B /;  target 
to  consist  of  a  compact  B  **  core  with  the  extra  neutron  loosely 
bound  to  it.  Y/hen  this  system  collides  with  the  deuteron,  it  is 
quite  conceivable  that  the  B  ,a  core  and  deuteron  should  merge, 
forming  a  C  /z* residual  nucleus,  while  the  loose  outer  neutron 
from  B  11  ,  only  distantly  involved  in  the  collision,  is  simply 

shaken  off  and  leaves  the  scene  of  action  with  the  momentum 

/ 

it  had  at  that  instant  (Figure  3). 
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This  mode  of  reaction  is  called  heavy-particle  stripping, 

9 

by  analogy  with  deuteron  stripping:  in  this  case  it  is  the 
heavy  particle,  (i.e.  the  B 11  ),  rather  than  the  deuteron,  which 
is  "broken  up  into  two  parts.  One  of  these  parts,  the  B  l°  core, 
hits  and  amalgamates  with  the  other  party  to  the  collision, 
namely  the  deuteron;  the  other  part  is  the  outgoing  neutron. 

For  comparison,  in  the  deuteron  stripping  mode,  it  is  the 
deuteron  which  is  "broken  up  into  two  parts:  one  of  these  parts, 
the  proton,  hits  and  amalgamates  v/ith  the  other  party  to  the 

s' 

collision,  namely  the  B ^  ;  the  other  part  is  the  outgoing 
neutron. 

The  idea  therefore  is  to  describe  heavy-particle  stripping 
by  the  usual  deuteron  stripping  theory,  simply  replacing  the 

deuteron  "by  the  B  ll  system,  its  proton  "by  the  postulated  3/0 

* 

core  of  the  B^  ,  and  the  31'  target  "by  the  deuteron. 

The  qualitative  reason  why  this  heavy-particle  stripping 
mode  is  suspected  of  being  responsible  for  the  anomalous 
angular  distributions  is  simply  this:  if  the  outgoing  neutron 
comes  from  the  B^  target,  and  its  final  momentum  is  compounded 
of  its  share  of  the  momentum  of  the  whole  target  plus  some 
contribution  from  its  internal  motion  in  the  B  “  ,  backward 
momenta  (  0^180*)  should  be  strongly  favored,  because  the 
target's  momentum  in  the  center-of-mass  system  before  collision 
is  in  the  direction  9  -  180°.  Just  as  deuteron  stripping 
implied  a  strong  forward  peak  in  the  neutron  angular  distribution 
because  the  neutron  was  originally  part  of  the  forward -moving 
deuteron,  heavy-particle  stripping  implies  a  strong  backward 
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peak  in  the  neutron  angular  distribution  because  the  neutron 
was  originally  part  of  the  backward-moving  heavy-particle  B11  . 
How  since  both  a  forward  and  a  backward  peak  occur  in  the 
anomalous  experimental  results,  it  ha3  been  concluded  (16)  that 
they  may  be  due  to  the  simultaneous  occurrence  of  both  modes  of 
reaction. 

The  idea  of  reactions  proceeding  by  more  than  one  mode  is 
of  course  not  special  to  this  reaction.  Because  of  the 
fundamental  indistinguis liability  of  elementary  particles  in 
quantum  mechanics,  it  is  never  certain,  in  atomic  or  nuclear 
reactions  involving  more  than  one  particle  of  the  sane  species, 
v/hether  the  observed  reaction  product  was  originally  part  of 
the  projectile  or  of  the  target.  The  standard  procedure  to  take 
this  into  account  is  to  antisymmetrize  all  wave  functions  vdth 
respect  to  the  coordinates  of  the  identical  particles. 
Observable  quantities  calculated  from  such  wave  functions. 

(e.g.  matrix  elements  or  reaction  cross-sections) ,  will 
therefore  contain  both  direct  and  exchange  terms.  The  direct 
terms  are  identical  with  the  ansv/er  that  would  have  been 

obtained  if  the  wave  functions  had  not  been  antisymmetrized. 

/ 

The  reason  why  the  antisymmetrization  can  often  be 

* 

neglected  in  practice  is  that  the  exchange  terms  turn  out  to  be 
small.  This  vd.ll  be  true  in  all  cases  where  the  observed 
reaction  product  "obviously"  has  a  much  greater  probability  of 
coming  from  one  of  the  reactants  than  from  the  other.  Por 
example,  ordinary  deuteron  stripping  depends  on  the  fact  that 
the  deuteron  is  a  very  loosely  bound  structure  a3  compared  to 
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other  nuclei.  Antisymmetrization  is  neglected  in  ordinary 
stripping  theory  "because  nucleons  within  the  target  nuclei  are 
normally  much  more  tightly  bound  than  in  the  deuteron,  and 
therefore  less  subject  to  the  type  of  collision  leading  to 
stripping,  mathematically  this  is  expressed  by  small  overlap 
between  the  neutron  wave  function  in  the  nucleus  and  that  in 
the  deuteron,  so  that  exchange  integrals  are  small. 

Thus  heavy -particle  stripping  is  always  present  in 
principle  but  normally  negligible  in  magnitude  compared  to 
deuteron  stripping.  The  question  is  whether  there  is  reason  to 
believe  that  in  this  particular  case  it  happens  to  be  comparable 
in  importance  with  deuteron  stripping. 

The  mathematical  manipulations  whereby  the  two  stripping 
modes  are  considered  in  the  evaluation  of  cross-sections  will 
be  reviewed  in  the  next  section.  This  will  be  followed  by 
sections  dealing  with  the  question  of  relative  magnitudes,  a.nd 
with  other  weak  points  in  the  two -mode  theory. 


G-  formal  Combination  of  the  two  modes 

The  evaluation  of  the  stripping  matrix  element  in  two-mode 
stripping  follows  the  particular  version  of  stripping  theory 
which  has  been  presented  in  Chapter  1.  The  important  new 
element,  the  antisymmetrization,  is  introduced  as  follows. 
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In  the  spirit  of  the  Born  approximation,  matrix  elements 
of  the  form 


will  he  considered.  One  of  the  two  wave  functions,  initial  or 
final,  must  he  antisymmetrized,  (it  can  he  shown  that  it  is  not 
necessary  to  anti  symmetrize  hoth  wave  functions  if  the 
interaction  potential  is  symmetrical  with  respect  to  the 
identical  particles.) 

B  11  can  he  assumed  to  contain  two  neutrons  deeply  hound  in 
the  Is  shell,  and  four  outer  neutrons  in  a  p  3/2  level.  These 
four  outer  neutrons  are  the  only  ones  considered,  on  the 
grounds  that  the  s  neutrons  are  much  less  likely  to  become 
separated  from  the  nucleus  in  a  stripping  type  of  collision. 

The  initial  wave  function  is  therefore  anti symmetrized  with 

/ 

respect  to. the  interchange  of  any  pair  of  the  five  neutrons 
which  may  turn  out  to  he  the  outgoing  one,  namely  the  one 
originally  in  the  deuteron,  and  the  four  originally  in  the 
p  3/2  outer  shell  of  the  B/;  nucleus. 

The  various  terms  in  the  wave  function  representing  the 
initial  state  will  contain  a  factor  representing  the  neutron 
from  the  deuteron,  and  a  factor  representing  the  four  neutrons 
in  the  p-shell  of  the  B  11  .  Since  only  the  neutrons  will  enter 
the  discussion  for  no w,  these  terms  can  he  abbreviated  like 
this : 

i  1,2345 


• 
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The  i  means  initial,  the  labelling  number  written  before  the 
comma  denotes  the  neutron  in  the  deuteron,  and  the  four 
labelling  numbers  v/ritten  after  the  comma  denote  the  neutrons 
in  the  p-shell  of  the  target. 

ITow  the  wave  function  describing  four  neutrons  in  this 
shell  must  itself  be  antisymmetric  with  respect  to  interchange 
of  any  pair  of  the  four  neutrons,  and  since  it  appears  as  a 
factor  in  the  various  terms,  thi3  property  greatly  simplifies 
the  antisymmetrization  that  must  be  performed  with  respect  to 
the  five  neutrons.  In  fact  an  antisymmetric  function  can  be 
written  in  just  five  terms: 


i  12345  /  c*  i  1,2545/- 


i  2,1345/ 


I i  3,2145 


-  I  i  4,2315  /  —  i  5,2341  > 


as  can  be  verified  by  taking  the  various  permutations  on  two 
particles . 

The  final  state  wave  function,  which  need  not  be  treated 


similarly,  can  be  written  f  1,2345  ,  where  f  stands  for 

final,  the  first  digit  represents  the  outgoing  neutron,  and  the 
four  digits  after  the  comma  represent  the  four  equivalent 

j,  A 

neutrons  left  in  the  p-shell  of  C  .It  is  of  course 
antisymmetric  under  interchange  of  any  pair  of  these  four 
indices  after  the  comma. 

The  interaction  potential  V  would  in  general  contain 
several  terms,  describing  interactions  between  neutron  and 
target,  proton  and  target,  and  neutron  and  proton,  in  the  case 
of  deuteron  stripping;  or  neutron  and  deuteron,  core  and 
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deuteron,  neutron  and  core,  in  the  case  of  heavy-particle 
stripping.  However,  (as  in  Section  B )  the  stripping 
approximations  consist  essentially  in  the  neglect  of  all 
interactions  but  the  one  between  the  tv/o  objects  which  actually 
hit  and  merge,  namely  Vp r  for  deuteron  stripping  vrhere  the 
proton  hits  and  is  absorbed  by  the  B  11  nucleus,  or  for 
heavy-particle  stripping  where  the  B/J  core  from  B  /l  hits  and 
coalesces  with  the  deuteron  (subscripts  p  for  proton,  T  for  the 

i 

B  ^  target,  d  for  deuteron,  and  C  for  the  core  of  B"  ). 

Thus  the  matrix  elements  of  these  interactions  are  to  be 
evaluated  between  a  five-term  wave  function  describing  the 
initial  state  and  a  single  term  wave  function  describing  the 
final  state,  in  order  to  give  a  transition  probability 

amplitude  T.  In  the  Johns  Hopkins  publications  the 

* 

simplification  of  the  matrix  element  by  means  of  the  stripping 
approximations  is  a  highly  formal  procedure,  involving 
expansion  of  the  five-term  matrix  element  in  terms  of  a  set  of 
initial  state  Green’s  functions,  and  neglect  of  all  but  the 
leading  terms  in  accordance  with  the  stripping  approximations 
relevant  in  each  case.  We  can  see  in  a  more  physical  way  the 
significance  of  these  manipulations.  We  have 


T  o<  <^f  1,2345  I  7  j  i  1,2345/'  -  f  1,2545  j  V  / 
-  /f  1,2345  |  V  J  i  3,2145  )>  -  <f  1,2345  jv  |l 

-  .  f  1,2345  V  i  5,254l\. 


i  2,1345y 
4 , 2315  -> 


How, 

interpret 


remembering  the  meaning  of  the  comma  notation,  we 
the  first  term  1,2345  jvj  i  1,2345  >:  here  the 
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neutron  which  was  originally  in  the  incident  deuteron  is  the 
one  labelled  1,  and  i3  the  one  that  goe3  out  and  is  observed: 
ordinary  deuteron  stripping  has  taken  place;  hence,  according 
to  the  stripping  approximations,  let  us  keep  only  Y ^ r  in  this 
case . 

But  in  the  second  term  /  f  1,2345  V  i  2,1345  ,  the 

\  »  I  / 

neutron  originally  in  the  incident  deuteron  is  the  one  labelled 
2,  and  it  is  not  the  same  as  the  outgoing  neutron:  it  goes  into 
the  while  the  neutron  labelled  1,  originally  in  the  3(/  , 

becomes  the  outgoing  particle.  Hence  we  say  that  heavy-particle 
stripping  has  taken  place  and  for  this  potential  V  we  keep  only 
;  -  ,  the  interaction  between  the  deuteron  (containing  the 
neutron  labelled  2)  and  the  B  1  core. 

Similarly  the  three  other  terms  \f  1,2345  J  V  j  i  3,2145.  , 

/t  1,2345  Jv  j  1  4,2315  )  and  f  1,2345  J  V  ,  i  5,2341 \ 
correspond  to  neutrons  3,  4,  and  5,  respectively,  in  the 
incoming  deuteron,  while  neutron  1,  originally  part  of  the 
target,  is  always  the  outgoing  one.  This  means  heavy-particle 
stripping;  hence,  by  the  stripping  approximations,  we  drop  all 
V  terms  but  Vj  ^  ;  ,  V  jfo)c  or  respectively. 

How  each  such  term  in  the  matrix  element  is  just  a  number,  the 
result  of  a  multiple  integration  and/or  summation:  the  labels 
1,  2,  ...  denote  different  variables  of  integration.  One  can 
therefore  relabel  at  will.  The  exchange  terms  are  now 
relabelled  so  that  in  each  case  the  neutron  from  the  deuteron 
gets  the  label  1,  and  the  outgoing  neutron  gets  the  label  2. 
this  end  the  operator  (^34^)  interchanges  the  labels  1  and  2 


To 
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/  |  /  \ 

in  the  first  exchange  term  v  i-2345  iva;  1 2>l34v>  > 

>  !  V 
giving  -  f  2,1345  ‘  Vy/;;£  ,  i  1,2345^  .  Sirci larly, 


1234 o \ 


f  .  X 


v 23145 


)  '  f  1,2345  .  Vv/i5r  i  3, 2145V  ">  <  f  2,3145  V  i  1,3245 


r.V 


(^2^42)<f  1,2345  V*  7  i  4,2315/ —>  f  2,4315  V..  i  1,4325 

V24315/  .  *  3  , 


(12345),^  lf2345(  , ->c  1  5,2341 


->  f  2,5341  V,  ,  ~  i  1, 


5342 


• 


Each  of  the  last  three  terms  in  the  form  obtained  above 


/ 


can  now  be  shown  to  be  equal  to  the  term  <  f 


O 

^  * 


1345 


V.  i 

l  »•/- 1 


1, 


2345 


by  means  of  interchanges  of  labels  within  the  p-shell  wave 
functions  only,  using  of  course  the  antisymmetry  property  of 
these  wave  functions.  It  can  readily  be  verified  that  whenever 
an  odd  number  of  interchanges  is  necessary  to  bring  the  initial 


wave  function  to  the  form  i  1,2345  ,  thU3  introducing  a 

*  / 

minus  sign,  another  minus  sign  cancelling  the  first  one  will 
similarly  be  produced  by  changing  the  final  wave  function  to 
the  form  i  f  2, 1345  . 

In  this  way  the  complete  matrix  element  for  two-mode 
stripping  reduces  to  the  interesting  form 


C-l 


T  o<  f 


1,2 


345. 


VT 

?T 


1,2345  /  -  4  f 

'  N 


2,1345  V 


i  1, 


2345 


> 
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This  completes  the  process  of  mixing  in  the  contributions  of 
the  two  modes.  The  first  term  is  given  by  pure  stripping  theory 
as  in  Chapter  1,  and  the  evaluation  of  the  second,  the  exchange 
or  heavy-particle  stripping  term,  follows  that  theory  as 
closely  as  possible,  with  appropriate  changes  of  course. 


. 
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If  one  is  interested  in  the  neutron  angular  distribution, 
the  cross-section  (Bl)  is  obtained  essentially  by  squaring  this 
matrix  element  T,  so  that  it  will  contain  a  deuteron  stripping 
term,  a  heavy-particle  stripping  term,  and  interference  terms. 
The  angular  correlation  between  the  neutron  and  gamma 
directions  is  obtained  by  inserting  the  matrix  element  T  of  G1 
into  D2. 


H-  Evaluation  of  the  matrix  elements 


In  two-mode  stripping  theory,  the  evaluation  of  the  matrix 
element  for  pure  deuteron  stripping,  <^f  1,234b  j  V^r  |i  1,234S^ 
is  performed  by  the  method  of  Chapter  1  but  presented 
differently,  and  further  assumptions  and  approximations  are 
made.  Since  the  physical  ideas  are  the  same  and  Edwards  1  text 
(6,7)  is  clear  on  this  point,  it  will  not  be  necessary  to 
introduce  the  complicated  notational  system  found  to  be 
necessary  to  carry  this  out  explicitly. 

The  procedure  is  to  carry  out  expansions  of  the  various 
factors  or  groups  of  factors  of  the  initial  state  wave  function 
in  terms  of  corresponding  factors  of  the  final  state  wave 
function,  using  the  physical  nature  of  the  process  as  a  guide 
in  choosing  the  proper  expansions.  Host  of  the  summations 
introduced  reduce  to  a  single  term  because  oi  the  expel imental 
conditions  (e.g.  selecting  only  states  with  a  given  outgoing 
neutron  momentum)  or  of  special  assumptions,  or  of  orthogomlity. 
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Thus  the  matrix  element  is  obtained  mainly  as  a  product  of 
coefficients  introduced  by  these  expansions,  such  as  various 
Clebsch-Gordan  coefficients,  factors  which  are  independent  of 
angle,  corresponding  to  the  factor  1^  in  the  cross-section 
formula  (Cl),  are  replaced  by  an  undetermined  proportionality 
constant. 

The  matrix  element  for  heavy-particle  stripping  is 
evaluated  in  the  same  way,  introducing  another  undetermined 
proportionality  constant .  According  to  Gl,  these  are  combined 
into  a  total  matrix  element  T  for  insertion  into  B1  or  D2, 
after  which  the  averaging  over  initial  and  summing  over  final 
quantum  numbers  must  be  performed. 

The  ratio  of  the  two  undetermined  proportionality 
constants  mentioned  above  is  not  evaluated  but  is  adjusted  to 
give  the  best  fit  to  the  experimental  data.  This  will  be 
discussed  below  (Section  j).  The  stripping  radius  for  each  of 
the  two  modes  is  also  kept  as  an  adjustable  parameter,  and  to 
obtain  the  good  fits  to  experimental  data  which  have  been 
claimed,  values  of  R  ranging  from  1.0  to  7.2  fermis  have  been 
quoted  in  various  JohnsIIopkins  publications.  These  three 
adjustable  parameters  do  not  influence  only  the  detailed 
structure  of  the  curves  obtained;  on  the  contrary,  their 
variation  within  ranges  of  values  as  reasonable  as  the  values 
actually  quoted  can  modify  the  appearance  of  the  curves  very 
considerably.  While  it  is  traditional  in  stripping  theory  to 
leave  one  continuous  variable,  the  radius,  as  an  adjustable 
parameter,  one  does  not  normally  accept  quite  that  much  leeway 
in  values  of  R* 
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The  other  ad  hoc  assumptions,  approximations,  and 
adjustable  parameters  introduced  will  now  be  listed,  to 
emphasize  the  very  considerable  arbitrariness  present  in  two¬ 
mode  stripping  theory.  It  should  be  borne  in  mind  that  all  this 
is  in  addition  to  the  assumptions  and  approximations  inherent 
in  stripping  theory  itself.  A3  discussed  elsewhere  in  this 
thesis,  the  criteria  for  their  applicability  in  this  particular 
case,  even  for  pure  deuteron  stripping,  are  hardly  satisfied  in 
a  convincing  manner;  their  application  to  heavy -imrticle 
stripping  is  untested  and  even  more  debatable. 

The  Bhatia  derivation  involved  a  sum  over  the  important 
angular  momentum  quantum  number  £ci  and  it  turned  out  (Section 
S)  that  only  the  values  -  l  and  3  are  allowed  by 
conservation  laws,  in  the  deuteron  stripping  treatment  of  this 
reaction.  Because  it  is  quite  reasonable,  in  shell-model 
terminology,  to  think  of  the  proton  as  going  into  the  available 
space  in  the  p  3/2  level  of  B  to  form  the  first  excited  C  , 
the  smaller  of  these  values  of  J(<  is  certainly  the  more 
important,  as  it  is  in  most  reactions.  Nevertheless  the  Born 
approximation  formula  includes  terms  corresponding  to  both 
values.  In  two-mode  stripping  theory  on  the  other  hand,  3uch  a 
summation  also  occurs  in  the  evaluation  of  tne  deuteron 
stripping  matrix  element,  but  it  is  arbitrarily  reduced  to  the 
single  term  J 1»  suppressing  the  ^=3  contribution 
entirely. 

This  approximation  seems  questionable.  Not  that  -he"  ^  is 
of  comparable  importance  to  /c  -  1  of  course  —  the  3he 11 -model 
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argument  is  sufficiently  decisive  on  this  point.  But  if  one  is 
going  to  take  into  account  the  heavy-particle  stripping  picture 
of  the  C  configuration  consisting  of  a  point  deuteron  with 
orbital  angular  momentum  0  or  2  around  a  3^  core,  the 
contribution  of  a  proton  with  (c~  3  may  not  be  negligible  in 
comparison  with  the  latter  picture,  particularly  since  the 
/c  =■  3  contribution  is  more  important  at  large  0  and  at  large 
radius  B  (as  can  be  seen  by  comparing  graphs  of  JJ  (K^)  j  *  (kR) 
for  various  R,  to  be  presented  in  Chapter  3);  the  deuteron 
stripping  radius  quoted  in  the  Johns  Hopkins  papers  (7)  is 
6.8  fermi3. 

The  same  sum  over  /C9  arbitrarily  reduced  to  one  term, 
occurs  also  in  the  evaluation  of  the  matrix  element  for  heavy- 
particle  stripping.  Here  the  angular  momentum  and  parity 
selection  rules  allow  / c~  4,  or  6.  This  refers  to  the 

orbital  angular  momentum  of  the  B  with  respect  to  the 
deuteron,  and  is  designated  as  /(C)  in  the  Johns  Hopkins 
papers.  The  fits  to  experiment  claimed  by  Ames  and  Owen  (17) 
and  quoted  by  Edwards  (7)  were  for  /(C)  —  0;  however  a  note 
added  in  proof  to  (7)  states  that  "Recent  experiments  of  Rask, 

Ame3 ,  and  Owen  [. .  .J  require  that  /(C)  =  2  rather  than  /(C)  -  1 

/ 

[(the  latter  is  a  misprint  for  /(C)  —  0)J  as  U3ed  in  the 
foregoing  It  is  obvious  that  this  quantity  i3  being 

used  as  yet  another  adjustable  parameter,  one  that  can  also 
produce  significant  changes  in  the  angular  distributions  and 
correlations . 
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Moreover  these  experimental  fits  require  an  assumed  spin 
of  1  for  the  3t0  core,  whereas  the  measured  spin  of  the  B >0 
nucleus  is  3.  If  such  an  assumption  is  justifiable  at  all,  it 
at  least  implies  that  this  quantity  is  also  being  used  as  an 
adjustable  parameter. 

The  evaluation  of  the  factor  ]T( Kb),  the  momentum 
transform  (C2  and  B7 )  of  the  deuteronfs  internal  wave  function, 
is  quite  approximate,  since  it  is  based  on  the  zero-range  wave 
function  only,  A  more  accurate  form  is  the  Hulthe'n  wave 
function 

KYr)-  jr  [  t"‘r -  cf'J 

> 

with  it  5,2  c<  (lloravesik,  20  ).  7T(K-b)  is  then  proportional  to 
£(K )J  X  instead  of  jU3t  (K^V  c(  x  )~2.The 

difference  is  unimportant  when  one  is  concerned  only  with  the 
position  of  the  principal  maximum,  but  becomes  significant  for 
a  study  of  the  whole  range  of  Q  from  0  to  180^,  especially  at 
2  Uev.  A  graph  of  this  quantity  Tf (K^)  in  arbitrary  units, 
plotted  as  a  function  of  k,  is  given  in  figure  3. 

liuch  more  uncertain  is  the  corresponding  factor  /T(K^) 
v/hich  appears  in  the  heavy-particle  stripping  matrix  element. 
This  is  the  momentum  transform  of  the  wave  function  for  the 
p-3hell  outer  neutrons  in  B  11  .  Owen  and  Uadansky  (12)  have 
calculated  this  assuming  a  square  well  potential,  and  this  is 
the  result  used  in  all  Johns  Hopkins  publications.  On  this 
point  it  is  possible  that  independent  experimental  evidence 
will  become  available,  from  measurements  of  (p,2p)  quasi-elastic 
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scattering  performed  with  high-energy  accelerators.  (Wattenberg 
21,  Wilcox  and  Eoyer  22).  These  experiments  yield  momentum 
distributions  for  protons  in  nuclei;  3ome  of  the  results  have 
been  interpreted  a3  showing  two  groups  of  protons,  namely  core 
protons  and  peripheral  protons.  Of  course  the  analysis  of 
these  experiments  involves  assumptions  (e*g.  the  impulse 
approximation),  and  the  deduction  of  neutron  momentum 
distributions  from  these  proton  momentum  distributions  will 
require  careful  justification;  but  it  is  possible  that  relevant 
information  may  come  out  of  this  experimental  method  in  the 
future,  after  a  larger  number  of  nuclei  are  studied  in  this  way. 

To  continue  the  listing  of  adjustable  parameters  in  two¬ 
mode  stripping  theory,  it  must  be  pointed  out  that  in  the 
evaluation  of  the  matrix  element  for  each  mode  one  has  to 
choose  between  the  various  coupling  schemes  possible  for  the 
angular  momenta  involved  in  the  reaction.  Y/hat  happens  is  that 
if  the  wrong  coupling  scheme  is  used,  expansions  are  made  in 
terms  of  a  set  of  eigenfunctions  of  a  quantity  which  is  not  a 
good  quantum  number,  and  hence  the  expansions  are  inconvenient. 

The  final  expressions  obtained  are  sensitive  to  the 
particular  coupling  scheme  chosen.  These  coupling  schemes,  even 
if  two  or  three  alternative  ones  are  considered,  are  themselves 
extreme  idealizations  of  the  real  physical  situation.  This  is 
recognized  by  Edwards  (7)  who  claims  without  proof  that  in  the 
general  case  of  a  mixed  coupling  scheme,  the  angular 
distributions  and  neutron-gamma  correlations  must  be  fitted  by 
linear  combinations  of  the  expressions  for  each  coupling  scheme, 
and  suggests  that  a  measure  of  the  amount  of  mixing  of  the 
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coupling  schemes  can  he  obtained  by  treating  the  coefficients 
in  these  linear  combinations  as  adjustable  parameters  to  fit 
the  experimental  results.  In  view  of  the  large  number  of  other 
adjustable  parameters  present  in  the  analysis,  this  suggestion 
must  not  be  taken  too  seriously. 

The  measured  angular  correlations  (lleilson  et  al.,  9)  can 
not  be  fitted  by  this  theory,  at  least  for  any  of  the  pure 
coupling  schemes. 


J-  Relative  laanitudcs 

9 

The  most  important  weakness  of  two-mode  stripping  theory 
is  that  the  parameter  expressing  the  relative  importance  of  the 
two  modes  is  left  open  for  adjustment  to  experimental  data. 

This  is  objectionable  because  variation  of  this  parameter 
within  a  range  of  values  which,  in  the  absence  of  any  published 
theoretical  estimate  Y/hafcever,  might  be  considered  acceptable, 
can  very  significantly  change  the  appearance  of  the  curves,  but 

r 

the  main  objection  to  thi3  procedure  is  the  following.  As 
discussed  in  Section  5*  of  this  chapter,  no  one  can  object  to 
heavy-particle  stripping  in  principle;  the  only  question,  the 
crucial  question,  is  whether  there  is  reason  to  believe  that  in 
this  particular  reaction  and  a  few  others  it  should  turn  out  to 
be  as  important  as  deuteron  stripping,  whereas  its  relative 
importance  is  iiormally  negligible.  By  leaving  this  relative 
importance  a3  an  open  parameter,  no  progre33  is  made  towards 
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answering  tills  question,  since  the  situation  is  otherwise  too 
indeterminate  —  involves  too  many  other  adjustable  parameters 
and  assumptions  —  to  say  that  the  relative  importance 
parameter  is  ndetermined,'by  this  best-fit  adjustment. 

Information  which  might  lead  to  estimates  of  the  relative 
importance  of  the  two  modes  is  very  sparse,  mainly  Decause  one 
does  not  normally  think  in  terms  of  point  deuterons  entering 
the  configurations  of  nuclear  states. 

Hadley  and  York  (23)  have  bombarded  a  C;i  target  with 
90  Hev  neutrons,  and  obtained  both  d  ~B/f  and  t  -f  3*  as 
reaction  products,  the  yield  of  tritons  being  at  most  10*.  of 
the  yield  of  deuterons.  This  can  be  explained  as  a  picx-up 
process  (Cnew  and  Goldberger,  24),  and  hence  a  careful  analysis 
might  3 how  to  what  extent  the  C/L  ground  state  can  be  thought 

/  Q 

of  as  consisting  of  a  deuteron  plus  B  core  rather  than  a 
proton  plus  B  ^  .  Unfortunately  we  are  concerned  with  an  excited 
state  of  C,Xf  rather  than  the  ground  3tate,  and  no  such  source 
of  information  i3  possible  for  excited  states. 

The  description  of  the  first  excited  state  of  C  tx  in  terms 
of  single  proton  excitation  from  the  p  3/2  to  the  p  £  level  is 
acceptable  but  far  from  unique.  The  collective  model  also 
accounts  for  ground  states  of  0"**  and  first  excited  states  ox  2 
for  even-even  nuclei.  Even  an  alpha -particle  model  for  C  4 
satisfactorily  correlates  the  ground  state,  the  2**  first 

a. 

excited  state,  the  0  level  at  7.65  Hev  and  the  level  at 

9 

9.61  Hev  (Glas3gold  and  Galonsky  25,  as  discussed  in  Almqvist 
et  al . ,  26). 


■ 
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Holmgren  et  al.  (27)  have  studied  the  reaction  B/,(lIe*,d)C/*  * 
(first  excited  3tate),  and  report  some  disagreement  with  the 
predictions  of  forward  stripping  theory.  They  state  that  the 
nature  of  the  additional  process  involved  can  not  he  predicted 
on  the  basis  of  their  present  results,  and  after  private 
communication  with  the  Johns  Hopkins  group  they  conjecture 
that  it  may  arise  from  exchange  stripping  or  compound  nucleus 
formation.  They  contend  that  “The  enhancement  of  this  process 
[heavy -particle  stripping]  relative  to  the  forward  stripping 
process  in  reactions  leading  to  the  first  excited  state  may 
possibly  be  due  to  the  ’cluster’  characteristics  (Wildernuth 
and  Kanellopoulos  £28  J)  of  the  low-lying  excited  states  of  C/a.  ” 

However  the  work  they  refer  to  (28)  does  not  support  this 
at  a.ll:  it  can  be  inferred  from  suggestions  made  in  (23)  that 
these  authors  (Wildermuth  and  Kanellopoulos )  expect  any  states 
involving  clusters  smaller  than  alpha  particles  to  be  of 
negligible  importance  in  C11.  Indeed  the  guiding  principle  of 
the  cluster  model  seem3  to  be  that  only  those  clusters  turn  out 
to  be  important  which  correspond  to  existing  nuclei  of 
exceptionally  tight  binding,  such »as  alpha  particles.  This  of 
course  argues  against  the  presence  of  a  deuteron  cluster  in  Ca7 
and  the  success  of  the  cluster  model,  really  an  alpha -particle 
model  in  this  case,  in  dealing  with  C,a  is  an  argument  against 
the  enhancement  of  heavy-particle  stripping. 

In  the  absence  of  anything  more  definite,  the  following 
classical^ argument,  suggested  by  Dr  Trainor  (29),  provides 
some  estimate  of  relative  importance.  Consider  the  break-up  of 
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either  the  deuteron  or  the  B ^  to  he  produced  by  the  impact  of 
their  collision.  The  work  done  by  this  impact  is  given  by  the 
product  of  the  change  in  momentum  by  the  change  in  velocity. 
Looking  at  the  problem  in  the  center-of-mass  system,  the  two 
colliding  systems  suffer  the  same  change  in  momentum;  but  the 
change  in  velocity  is  much  more  severe  for  the  deuteron  than 
for  the  B  ^  ,  (by  a  factor  of  11/2) ,  because  of  its  smaller 
mass.  Thus  by  times  as  much  of  the  work  done  by  the  impact  of 
the  collision  i3  available  to  the  deuteron,  for  breaking  bonds, 
as  to  the  B  11 . 

In  the  case  of  the  deuteron  this  energy  of  impact  acts 
entirely  on  the  single  bond,  but  in  the  case  of  B  ^  it  is 
divided  among  at  least  the  four  bonds  holding  the  neutrons  we 
are  interested  in;  hence  twenty-two  times  as  much  energy  is 
available,  on  the  average,  for  breaking  the  deuteron  bond  &3 
for  breaking  any  one  of  the  four  p-shell  neutrons’  bond  in  B^  • 

These  quantities  of  energy  must  be  compared  with  the 
binding  energy  in  each  case.  Thi3  introduces  another  factor  of 
5.2  in  favor  of  the  deuteron,  since  the  binding  energy  of  the 
last  neutron  in  B  rl  is  11.5  Lev  a3  compared  to  2.226  Lev  for 
the  deuteron.  Thus  at  2  Lev  bombarding  energy  for  example, 
kip.  Ay  z  h*k(£/m  is  about  2.38  Mev  for  the  deuteron,  compared 
to  a  binding  energy  of  2.226  Lev;  whereas  it  is  only  about 
.524  llev  for  B11  ,  or  at  most  .13  Lev  per  bond  on  the  average, 
compared  to  a  binding  energy  of  11.5  Lev'.  ±Ience  a  Break-up  ox 
B  11  leading  to  heavy-particle  stripping  is  much  less  probable 
than  a  break-up  of  the  deuteron  leading  to  deuteron  stripping. 
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In  addition  to  this  enormous  advantage  for  deuteron 
stripping  relative  to  heavy -particle  stripping,  there  must  also 
be  an  additional  advantage  coming  from  the  reduced  widths  of 
the  final  state,  since  B  11  +  p  is  certainly  a  better  picture  of 

I  1  *¥  IQ 

C  than  B  +  d,  even  if  it  is  difficult  to  specify  how  much 
better.  Admittedly  these  calculations  are  based  on  a  very 
inadequate  model  of  the  stripping  process,  but  at  least  they 
illustrate  why  one  is  reluctant  to  consider  heavy-particle 
stripping  to  be  as  important  as  deuteron  stripping. 

The  general  conclusion  that  3eem3  to  emerge  from  all  this 

is  that  no  valid  reason  whatever  has  been  put  forward  to  expect 

the  relative  importance  of  heavy-particle  stripping  relative  to 
deuteron  stripping  to  be  greater  in  this  reaction  and  a  few 
others  than  it  normally  is;  and  that  the  analysis  performed  on 
this  assumption  involves  too  many  approximations  and  adjustable 
parameters  for  the  fits  to  experiment  to  be  considered  as 
verifying  it. 

This  does  not  detract  from  the  value  of  the  work  reviewed 

in  shOY/ing  how  the  complicated  kinematics  may  be  handled  and 

the  two  modes  formally  combined.  The  analysis  may  be  of  value 
in  another  reaction,  or  under  experimental  conditions  such  that 
one  can  be  3ure  that  no  trouble  arises  from  Coulomb  or  nuclear 

effects  neglected  by  the  basic  stripping  approximations. 

/ 

q 

Unfortunately  the  Be  nucleus,  which,  because  of  the  very 
low  binding  energy  of  it3  last  neutron  (1.7  Lev),  would  3eem  to 
be  the  most  favorable  subject  for  a  heavy -particle  stripping 
study,  i3  likely  to  be  very  complicated  to  3tudy  because  it 
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geems  that  one  should  really  consider  it  as  two  alpha, -particle 

a 

clusters  plus  an  extra  neutron,  rather  than  a  Be  core  plus  a 

neutron.  loreover,  the  work  of  ITeilson  et  al.  (9)  on  this 
nucleus  shows  angular  distributions  which,  for  all  hut  one  of 
the  states  of  B /0  ,  do  not  resemble  the  curves  predicted  by 
heavy-particle  stripping  theory. 

In  high-energy  bombardment  of  deuterons  by  deuterons,  of 
course,  two  modes  must  be  considered. 

finally,  one  must  remember  in  this  type  of  work  that  the 
interference  term  may  be  quite  large  even  though  the  pure 
heavy-particle  stripping  term  is  negligible. 


. 
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Chapter  5  ALTEHhAlIVB  EEPIAEAII0H3 


IZ-  qualitative  Discussion 


The  special  features  of  the  neutron  angular  distributions 
observed  experimentally  for  By/  (d,n)C/j<  (4.43  Eev  state)  and 
other  similar  reactions  have  been  described  in  Section  E,  and 
some  shortcomings  of  the  heavy -particle  stripping  proposal  as  a 
detailed  explanation  of  them  have  been  indicated  in  Chapter  2. 

In  the  present  work,  a  study  has  been  made  of  the  possibility 
of  accounting  for  these  special  features?  in  particular  the 
large  cross-sections  for  0  j>  90° ?  without  introducing  the 
heavy -particle  stripping  hypothesis. 

It  is  not  contended  that  heavy-particle  stripping  doe3  not 
take  place,  for  in  principle  of  course  such  effects  must  be 
considered.  However  the  work  of  the  Johns  Hopkins  group 
(Chapter  2)  which  attempts  to  explain  the  backward  distributions 
as  a  heavy-particle  stripping  effect  suggests  that  heavy- 
particle  stripping  contributes  equally  with  deuteron  stripping. 
In  the  present  chapter,  the  point  of  viev/  is  adopted  that 
heavy-part icle  stripping  contributes  only  in  a  minor  way,  and 
that  one  must  seek  elsewhere  for  an  explanation  of  the  backward 
distributions . 

One  suggested  alternative  explanation  can  be  summarily 
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dealt  with.  Because  some  results  for  3i  (d,p)Si  and 
A1  27(d,p)Al 18 ,  also  showing  abnormally  high  cross-section  at 
backward  angles,  3eemed  well  fitted  by  stripping  curves  with 
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the  factor  77"( Kf>)  replaced  "by  a.  constant,  Bromley  et  al.  (44) 
have  examined  a  number  of  deuteron  potentials,  and  demonstrated 
that  no  reasonable  potential  will  produce  such  a  constant  form 
factor. 

There  is  reason  to  think  that  Coulomb  effects,  not 
considered  at  all  in  the  basic  theory  of  deuteron  stripping, 
may  be  particularly  important  in  this  reaction.  Criteria 
designed  to  establish  whether  their  neglect  is  justified  in 
particular  cases  (Butler,  1)  are  not  satisfied  for  this 

reaction  at  the  bombarding  energies  used,  3ince  the  Coulomb 

/ 

barrier  a3  seen  by  the  approaching  deuteron  is  of  the  same 
order  of  magnitude  as  these  bombarding  energies.  The  reason  one 
tries  to  fit  results  on  the  basis  of  stripping  theory  at  all  is 
that  in  other  cases  it  often  seems  to  work  reasonably  well  even 
when  the  criteria  are  not  satisfied.  The  procedure  of  keeping 
the  stripping  radius  as  an  adjustable  parameter  to  obtain  the 
best  fit  to  experiment  seems  to  cover  a  multitude  of  sins.  But 
as  the  particular  results  under  study  turn  out  to  be  very 
different  from  the  usual  stripping  distributions,  it  is  natural 
to  a3k  whether  one  can  see  in  a  simple  way  if  Coulomb  effects 
might  modify  the  distributions  in  the  required  way. 

qualitatively,  this  is  confirmed  by  descriptive  arguments, 
which  vary  depending  on  how  one  chooses  to  look  at  the  proce33. 
Bor  a  (d,n)  reaction  the  standard  qualitative  description  (e.g. 
Butler  (l),  Chapter  l)  considers  separately  the  effect  of 
Coulomb  forces  on  the  trajectory  of  the  approaching  deuteron, 
and  on  the  wave  function  describing  the  proton  after  capture  by 
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the  nucleus.  The  latter  effect  can  he  seen  to  lead  to  a 
considerable  reduction  in  the  total  cross-section  hut  has 
little  effect  if  any  on  the  angular  distribution.  It  can  he 
thought  of  a3  acting  only  on  the  capture  probability  of  Cl, 
which  is  independent  of  angle,  and  is  therefore  of  no  interest 
for  present  purposes.  The  Coulomb  effect  on  the  deuteron 
trajectory,  on  the  other  hand,  is  interesting,  and  even  a 
purely  classical  description  of  it  brings  out  the  essential 
features  of  the  modifications  thus  introduced. 

Classically  a  deuteron  approaching  with  kinetic  energy  E 
a  center  of  repulsive  Coulomb  force  Ze^/r2'  at  an  impact 
parameter  b,  suffers  an  angular  deflection  of  6  away  from  the 
incident  direction,  where 


tan 


If  a  stripping  process  occurs  which  would  otherwise  have 
produced  an  outgoing  neutron  at  an  angle  0 ,  the  result  roughly 
is  an  outgoing  neutron  at  an  angle  0  ~  0  f  ^  ,  since  the 
Coulomb  field  has  already  led  to  a  deflection  ^  . 

A  good  idea  of  the  resulting  neutron  angular  distribution 
can  therefore  be  obtained  by  Just  replacing  0  by  0  in  the 
uncorrected  angular  distribution,  with  b  of  the  order  of  the 
stripping  radius  R  or  somewhat  greater.  Thus  peaks  and  minima 
of  the  angular  distribution  are  shifted  to  larger  scattering 
angles  by  a  constant  amount  when  Coulomb  forces  are  taken  into 
account  in  this  crude  manner,  the  amount  decreasing  as  the 
bombarding  energy  increases  relative  to  the  Coulomb  barrier 
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height.  Moreover,  since  many  different  impact  parameters 
simultaneously  contribute,  it  is  evident  that  the  sharp 
features  of  the  stripping  pattern  are  smeared  out  or 
attenuated  to  some  extent. 

How  the  observed  anomalous  angular  distributions  for 
3  (d,n)C  ^  might  be  described  as  an  extreme  case  of  3uch  a 

shift  in  pattern  and  attenuation  of  features.  Of  course  the 
classical  argument  above  can  only  indicate  the  necessity  of  a 
proper  calculation  of  the  Coulomb  effect.  This  turns  out  to  be 
extremely  difficult  to  do. 


L-  Coulomb  functions 

This  section  introduces  the  mathematical  functions  which 
have  been  found  necessary  both  in  the  standard  treatments  of 
Coulomb  effects  in  stripping  (Cection  M)  and  in  the  different 
treatment  proposed  here  (Section  H) . 

The  approach  of  a  particle  of  reduced  mass  m  and  charge 
-he,  with  momentum  k  =  mv  with  respect  to  a  scattering  center  of 
charge  Ze,  is  governed  by  the  Schrfldinger  equation 


A  solution  of  this  equation  which  describes  the  scattering  in  a 
physically  acceptable  way  can  be  found  (llott  and  La33ey,  30)  by 


making  the  substitution 

that  F  ^F(r-z).  The  Frobenius  method  of  solving  the 


fact 
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differential  equation  then  leads  to 


L2 


F  (-'?■,  I)  iho-cos&i) 


where 


Ze 


t  V 

and  g  (x)  is  the  gamma  functic 


7fKAr 


ion.  The  function 


a  v  ,  a  (a+t)  z 

b  b  (b+t)  2l 


is  called  the  confluent  hypergeometric  function  (or  aeries)  by 
many  authors  (jahrikre  and  Smde ,  31;  Bloch  et  al.,  32),  and  the 
ICummer  function  by  others  (Bagnus  and  Oberhettinger ,  35). 

It  is  more  useful  for  many  purposes  to  solve  LI  by 
separation  of  variables  in  spherical  polar  coordinates .  This 
procedure  leads  to  a  radial  equation 


A  solution  of  thi3  can  be  written  in  terms  of  confluent 
hypergeometric  functions  different  from  those  used  above: 


iFt  (‘7+  F+l;  -tiki) 


These  solutions  L^(r)  occur  as  coefficients  in  the  very 
important  spherical  harmonics  expansion  of  


established  by  Gordon  (34): 


L4 


where  (5^  =  arg  p  (  /+  1  r  i  7  )  • 
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One  prefers,  nowadays,  to  transform  the  radial  equation  L3 
by  introducing  the  new  variables  kr  and  u  ((’)  =  pup)- 

Thi3  gives 


cTu 

<*r 


Two  independent  solutions  of  this  last  form  of  the  equation  are 
known  as  the  Coulomb  functions.  They  are:  the  regular  function 
a  real  function  which  vanishes  at  f  —  0  and  i3 
normalized  30  that  its  asymptotic  behaviour  for  f7  is 


sin 


frr 

Z 


rA  v-  p) 


1 

J 


and  the  irregular  solution  Gj>[  p)  which  doe3  not  vanish  at 
—  0,  and  has  the  corresponding  asymptotic  form 

G,(f)  =  If- *£-7 *'<>■()  +  <$]  • 

This  normalization  implies  that  the  Y/ronskian 

<iZ  f  d  6j>  / 

df  "  £  ~df  =  ' ' 

Ty(  (>)  and  (  p)  are  of  course  functions  of  the  variable  7 
which  measures  the  strength  of  the  Coulomb  interaction. 


k-  Coulomb  Effects  in  Stripping 

The  introduction  of  the  effects  of  Coulomb  forces  into  the 
standard  theory  of  stripping  angular  distributions  ha3  been 
studied  by  many  authors  (e*g.  34  to  39).  The  following  features 


are  typical  of  such  treatments. 
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Calculations  are  performed  for  the  (d,p)  reaction,  a  brief 
footnote  usually  claiming  without  discussion  that  they  can  be 
applied,  mutatis  mutandis,  to  the  (d,n)  reaction  also. 

The  method  is  simple  in  principles  the  plane  waves 
describing  the  approaching  deuteron  (considered  a3  a  point 
particle)  and  the  outgoing  proton  are  replaced  by  the 


(see 


appropriate  Coulomb-distorted  \?ave  functions 


previous  section),  and  the  analysis  is  carried  through  along 
the  general  lines  of  the  non-Coulomb  treatment,  except  that  the 
simplicity  and  significance  of  the  derivation  as  presented  in 
Chapter  1  are  lost  in  the  process. 


These  functions 


instance  as  contour  integrals  in  the  complex  plane,  whose 
arguments  are  expanded  in  various  ways,  leading  to  a  formidable 
amount  of  numerical  analysis  and  computation,  much  of  which 
must  be  performed  anew  for  each  particular  reaction  studied. 

These  calculations  have  been  performed  for  several 
particular  cases,  and  the  results  obtained  generally  do  confirm 
the  three  trends  brought  out  qualitatively  in  Section  K,  i.e. 
reduction  in  absolute  cros3-section,  shift  of  the  recognizable 
stripping  pattern  towards  larger  $,  attenuation  of  the 
sharpness  of  its  features;  and  they  generally  fit  experimental 
results  somewhat  better  than  the  uncorrected  formulas. 

These  calculations  treat  the  deuteron  as  a  point  charge, 
rather  than  recognizing  the  fact  that  the  Coulomb  field  really 
acts  on  the  proton.  Some  treatments  attempt  to  alleviate  this 


in  an  uncertain  way  by  using  a  wave  function  which  represents  a 


. 
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particle  having  the  deuteron’s  mass  and  charge  "but  the  position 
coordinates  of  the  proton. 

A  full  treatment  of  the  stripping  problem  would  also 
include  corrections  for  the  complex  nuclear  interactions 
involved,  which  are  largely  neglected  "by  the  "basic  stripping 
assumptions.  This  can  not  "be  done  adequately  3ince  nothing  but 
a  few  generalities  are  known  unequivocally  about  these 
interactions,  as  well  as  about  detailed, nuclear  structure. 

A  few  calculations  have  been  made  which  try  to  take  the  nuclear 
interactions  into  account  by  introducing  simple  models,  for 
particular  cases,  in  combination  with  Coulomb  effect 
calculations  mentioned  above;  the  same  comments  apply.  The 
problem  reduces  in  principle  to  a  determination  of  the  amount 
of  irregular  Coulomb  function  to  be  properly  mixed  in 

with  the  regular  Ty (  p ) ;  in  practice  this  tends  to  complicate 
the  already  severe  computational  problem  by  orders  of  magnitude 
considering  the  various  possible  models  to  be  considered  and 
the  additional  unknown  parameters  involved.  The  trend  of 
results,  for  those  cases  which  have  been  considered,  i3  to 
cancel  out  part  of  the  pure  Coulomb  modifications  outlined 
above.  For  practical  purposes  it  seems  that  a  good  case  can  be 
made  for  considering  Coulomb  corrections  while  still  neglecting 
nuclear  corrections,  unless  one  is  very  close  to  a  compound 
nucleus  resonance  level.  (Resonances  do  occur  in  reactions 
which  show  good  stripping  angular  distributions.) 
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II-  Proposal  for  (d,n)  Reactions 

,  j?rom  the  point  of  view  of  trying  to  explain  the 
B  ^  (d,n)C  /a  ^"anomalous  angular  distributions  a3  a  special 
Coulomb  effect,  the  work  mentioned  in  the  last  section  is  not 
very  helpful.  The  computations  are  too  complicated.  It  was 
designed  to  deal  with  cases  where  the  characteristic  stripping 
pattern  is  somewhat  modified,  not  essentially  transformed.  And 
it  was  designed  to  deal  with  (d5p)  reactions  rather  than  (d,n). 

The  last  point  is  particularly  significant.  It  may  be 
adequate  to  consider  Coulomb  effects  on  a  point  deuteron  when 
one  is  discussing  (d,p)  stripping,  in  which  the  proton  never 
comes  close  enough  to  the  target  nucleus  to  feel  the  full  force 
of  the  Coulomb  interaction,  although  there  i3  some  current 
uneasiness  about  this.  However,  there  i3  much  less  justification 
for  this  approach  in  (d,n)  reactions,  particularly  if  one  is 
looking  specifically  for  significant  modifications  of  the 
stripping  angular  distributions  due  to  Coulomb  effects  on  the 
close  approach  of  the  charged  proton. 

This  should  be  made  more  forceful  by  another  version  of 
the  qualitative  argument  which  leads  U3  to  expect  that  a  strong 

t  As  evidenced  for  instance  by  thi3  statement  from  the  latest 
major  work  on  the  subject  (Satchler  and  Tobocman,  39):  "  it  is 
not  certain  even  that  merely  distorting  the  center-of- 

mass  motion  of  the  deuterons  is  an  adequate  representation  of 
their  wave  function  close  to  the  nucleus." 


■ 
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Coulomb  effect  can  favor  large  angle  stripping  in  a  (d,n) 

/ 

reaction.  Consider  the  problem  in  terms  of  barrier  penetration. 
Barrier  penetration  by  the  deuteron  is  not  a  very  helpful 
concept  in  this  context  because  there  is  no  unique  ra.dius  to 
which  the  deuteron  must  penetrate  for  stripping  to  take  place. 
But  barrier  penetration  by  the  proton  is  more  informative , 
since  the  proton  must  go  in  far  enough  to  be  stripped  off  by 
nuclear  forces.  Bow  the  probability  of  penetration  of  a  Coulomb 
barrier  is  a  very  sensitive  increasing  function  of  the  momentum 
of  approach.  The  presence  of  the  barrier  thus  favors  high 
momentum  of  approach  k  and  hence  large  angle  0 .  It  is  obvious 
that  to  study  3uch  an  effect  quantitatively  one  must  consider 
Coulomb  effects  on  the  proton’s  approach  to  the  nucleus  with 
momentum  k  rather  than  on  the  whole  deuteron. 

We  are  dealing  here  with  a.  quantum  mechanical  three-body 
problem  involving  proton,  neutron,  and  target,  which  must  be 
reduced  to  a  more  amenable  problem  by  one  trick  or  another.  The 
work  appearing  in  the  literature  to  date  has  reduced  it  first 
to  a  tv/o-body  problem  involving  point -deuteron  and  target.  3uch 
approaches  do  not  do  justice  to  the  physical  process  described 
in  the  previous  paragraph.  An  alternative  approach  to  the 
reduction  of  the  three-body  problem  is  considered  here.  The 
point  is  that  in  the  case  under  study,  the  deuteron  bombarding 
energy,  the  binding  energy  of  the  deuteron,  and  the  Coulomb  and 
centrifugal  potentials  of  consequence  between  proton  and  target 
are  all  of  the  some  order  of  magnitude,  about  1-  to  5  Bev,  so 
that  alternative  reductions  may  be  equally  justifiable. 
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Consider  the  three-factor  formula  Cl 

<r  =  /J 

Sc  ^ 

discussed  extensively  in  Section  C.  Ignoring  entirely  Coulomb 
effects  on  the  deuteron  as  such,  let  us  try  simply  to  replace 
the  factor  Ly  [k),  describing  the  approach  of  the  proton 

towards  the  target  a3  a  plane  wave  with  momentum  k,  by  a 

/ 

corresponding  factor  L^?(!c)  describing  as  nearly  a3  possible 
the  same  situation  with  Coulomb  forces  included. 

The  factor  L^(k)  is  the  significant  one,  which  determines 
the  characteristic  features  of  the  usual  stripping.  It  seems 
reasonable  that  modification  of  the  factor  L  ✓  (k )  should 
account  for  the  major  Coulomb  effects  on  the  close  approach  of 
the  proton  leading  to  its  absorption  into  the  target,  whereas 
the  arguments  advanced  to  justify  the  standard  Coulomb 
treatment  of  (d,p)  reactions  (Section  ll)  still  have  3ome 
validity  for  description  of  the  distant  approach  of  the  whole 
deuteron.  These  arguments  hold  that  because  the  Coulomb  field 
does  not  change  too  drastically  during  the  time  corresponding 
to  the  internal  period  of  the  deuteron,  the  deuteron*s  internal 
motion,  and  hence  Tf  (IC^  ) ,  is  not  too  affected  by  them.  Coulomb 
effects  on  the  factor  iy  do  not  affect  the  angular 
distribution,  and  can  thu3  be  left  out. 

The  appropriate  form  of  L^(k)  has  been  deduced  by  antilogy 
with  the  simple  Born  approximation  formula  developed  and 
discussed  extensively  in  Chapter  1.  There  L^.(k)  arises  as  the 
modulus  squared  of  the  coefficient  in  the  spherical  harmonics 
expansion  B8  of  a  plane  wave  representing  the  approach  of 
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proton  to  target  with  relative  momentum  k.  Retaining  as  much  as 
possible  of  this  picture  of  the  break-up  of  the  deuteron 
producing  such  a  plane  wave,  we  consider  what  is  the  nearest 
thing  to  a  plane  wave  possible  in  the  presence  of  Coulomb 
forces.  This  is  the  function  (L2) 

%  ft  $)  =  a*6  <•*&). 


For  the  particular  case  of  vanishing  Coulomb  forces, 
y  -  0,  the  differential  equation  LI  satisfied  by  ^  (r ,  0) 
reduces  to  the  free-particle  Schrddinger  equation 

\7X?  r  k'H'  =  o 

ike  cos  G 

of  which  the  plane  v/ave  €  is  a  solution,  and  each  of 

the  factors  0  ro^7)  and  /  F,  reduces  to  unity, 

,  _  ,  _  n  (k'e  COS  & 

leaving  only  the  plane  wave  0 

iloreover  the  Gordon  expansion  L4  of  ^(r,  0),  in  modern 

notation 


corresponds  precisely  to  the  expansion  138  of  a  plane  wave  in 
spherical  harmonics.  The  (  C)/^  are  solutions  of  the  rokiial 
Coulomb  equation 


which  reduces,  for  ~  0,  to  the  equation 


A 
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t*-  dr 


-  [k' - 'COM)  L  -o 
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of  which  the  spherical  Bessel  functions  i/e  )  appearing  in  the 
uncorrected  stripping  formula  are  solutions.  It  follows  that 
the  required  factor  is 

Lj,  W  =  tr  (zA  1)  /// 

In  addition  to  ( f>)/r  ,  equation  in  also  has  other 
solutions  G  r^fvf  ,  irregular  at  the  origin,  just  as  equation 
1T2  ha3  solutions  n^(  (9 )  in  addition  to  SA/»  .  These  are 
excluded  from  the  study  of  many  idealized  physical  problems 
because  they  diverge  at  the  origin  and  are  thus  not  physically 
acceptable.  However  in  problems  where  the  range  of 
applicability  of  this  differential  equation  does  not  include 
the  origin,  this  objection  does  not  arise  and  the  general 
solution  includes  some  contribution  from  the  irregular  solution 
as  well,  the  relative  amount  depending  on  the  boundary 
conditions.  The  problem  of  a  particle  approaching  a  nucleus  is 
in  fact  such  a  problem,  since  the  Schrddinger  equation  for  the 
interior  of  the  nucleus,  which  includes  the  origin  of 
coordinates,  is  different  from  the  simple  Coulomb  equation,  and 
a  fully  rigorous  solution  (along  the  lines  of  the  Butler 
approach  if  it  did  not  involve  subsequent  approximations)  would 
certainly  involve  the  irregular  solutions.  But  the  Born 
approximation  treatment,  neglecting  the  effect  on  the  plane 
wave  of  the  presence  of  the  nucleus,  gives  good  results  in  the 
uncorrected  stripping  case  without  involving  the  irregular 
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spherical  ITeumann  functions  n^(  ^  ) ;  it  therefore  3eems 
reasonable  to  U3e  this  L g  (k)  involving  only  the  regular 
The  proposed  formula,  then,  is 


with 


2. 


replacing  C3 


1  /c 


113 

1T4 


P-  A  Uethod  of  Hxamining  Reaction  Data 

In  the  course  of  this  work  a  method  of  examining 
experimental  angular  distributions  has  been  used  which  can  be 
helpful  and  timesaving  in  deciding  whether  certain  data  exclude 
a  pure  deuteron  stripping  explanation  or  some  simple 
modification  of  it  such  as  the  one  suggested  above.  The  method 
does  not  seem  to  have  been  pointed  out  explicitly  before, 
although  the  approximations  on  which  it  is  based  are  generally 
known  and  used  in  numerical  calculations. 

Again  we  use  the  three-factor  cross-section  formula.  Cl, 

<r  -  TT(K)  2  (V  £ 

and  neglect  which  is  independent  of  angle.  The  most 
important  factor,  L/?  (k),  depends  on  angle  and  on  bombarding 
energy  only  insofar  as  it  depends  on  lc,  which  is  given  by 


■ 
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^  '  4  +  k)  -  £ 

As  to  the  factor  77*(K^)  ,  its  dependence  on  G  and  on 
energy  is,  strictly  speaking,  through  the  quantity  v/hich  is 
determined  "by 

~  k  +  i  /cy  —  k  ^dcos  ® 


But  Ki  is  related  by  conservation  laws  to  the  momentum  k,  and 
it  turns  out  that  this  relation  is  given  to  a  good 
approximation  by  a  formula 


+-  ^ 


XM/x 


av  xx)  . 


which  does  not  otherwise  involve  bombarding  energy  or  angle. 
I.e.,  for  a  given  reaction  to  a  given  level  of  the  final 
nucleus,  all  the  quantities  appearing  in  this  approximate 
formula  except  and  k  themselves  are  independent  of  angle  and 
of  bombarding  energy.  (  o(  is  the  internal  wave  number  related 


to  the  binding  energy  of  the  deuteron 


~h  xo< 


,  with  m  =- 


m. 


and  7C  is  similarly  related  to  the  binding  energy  jL. 3T  of  the 

Zmc 

captured  nucleon  c  in  the  final  nucleus.)  Thus  the  factor 
IT (K#  )  also,  is  expressed  for  purposes  of  computation  as  a 
function  of  k,  which  doe3  not  depend  on  angle  except  insofar  as 
it  depends  on  k.  This  is  therefore  true  of  the  whole  relative 
angular  distribution  expression. 


Derived  from  an  equation  amounting  to  conservation  of  energy, 
by  setting  masses  of  nucleons  and  nuclei  (in  atomic  mas3  units) 
equal  to  the  nearest  integer. 
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Suppose  now  that  experimental  data  for  the  same  reaction 
(same  level)  at  different  bombarding  energies  is  available.  If 
the  stripping  formula  is  to  describe  it,  the  variation  of  G~ 
with  0  at  different  energies  should  be  such  that  the  cross- 
section  as  a  function  of  k  be  unique.  The  cross-section  at  one 
combination  of  energy-value  and  angle-value  giving  a  certain  k 
should  be  the  3ame  as  that  at  any  other  combination  of  energy- 
value  and  angle-value  giving  the  same  k.  Experimental  curves  of 
cross-section  as  a  function  of  0  for  different  energies  should 
therefore  coalesce  to  a  single  curve  when  the  points  are 
replotted  as  a  function  of  k  instead  of  0  .  Or  at  least  their 
deviations  therefrom  should  be  random,  not  systematic. 

Before  applying  this  test  the  curve  for  each  energy  may  be 
renormalized  ad  hoc  since  the  cros3-section  measurements  are 
only  relative  in  any  case.  The  variation  of  with  bombarding 
energy  is  therefore  irrelevant.  Horeover  an  arbitrarily 
adjustable  amount  of  isotropic  background  blamed  on  processes 
other  than  stripping  may  be  subtracted  from  the  experimental 

9 

curves. 

If  the  data  are  such  that  they  can  be  fitted  by  the  pure 
stripping  formula  Cl  (either  Born  approximation  or  Butler),  for 
any  R,  and  for  any  combination  of  values  of  3^  ,  (i.e.  even  for 
the  unexpected  case  of  3E?<»  3^  where  1  );  or  by  the 

modified  formula  113  embodying  Coulomb  effects  as  suggested  in 
the  preceding  section;  or  by  a  similar  formula  designed,  to  cope 
with  other  modifications  within  the  framework  of  the  very  general 
formula 
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<r  - 


irofv  ;>  Q  m  ft. 


a 


as  long  as  Ly  (k)  is  a  function  of  6  only  through  k;  then  a 
unique  curve  for  <7"(k)  will  be  obtained,  and  the  task  of 
looking  for  a  theoretical  expression  to  predict  it  should  he 
facilitated  by  that  uniqueness.  The  possibilities  mentioned 
imply  such  a  great  variety  of  possible  shapes,  considering  the 
great  latitude  currently  fashionable  in  the  choice  of  R,  that 
merely  inspecting  curves  of  <T  {  6)  Y/ould  not  have  helped  to 
decide  this. 

If  on  the  contrary  the  data  include  a  significant  non¬ 
isotropic  contribution  from  processes  other  than  deuteron 
stripping,  such  as  for  instance  heavy-particle  stripping,  it  is 
very  unlikely  that  they  should  exhibit  such  uniqueness.  One  has 
then  a  valuable  indication  that  such  processes  must  be  taken 
into  account. 

The  procedure  then,  is  first  to  calculate  k/  and  ly  for 
the  various  experimental  bombarding  energies,  from  the 
energetics  of  the  reaction,  then  to  tabulate  or  graph  k  as  a 
function  of  6  (from  Equation  11)  for  each  energy. 

It  7/ill  be  seen  from  PI  that  the  range  0  ^  9  £.  180° 

corresponds  to  greater  and  greater  ranges  of  k  as  the 
bombarding  energy  increases;  the  range  of  k  at  a  certain 
energy  includes  the  ranges  obtained  for  all  lov/er  energies. 
Using  larger  and  larger  bombarding  energies  allows  us  to  look 
at  a  greater  and  greater  portion  of  the  curve  of  (T  as  a 
function  of  k. 
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One  3hould  then  replot  each  curve  C  ( 0  )  as  ^“(k)  and 
renormalize  so  that  the  curves  for  different  energies  will  all 
agree  at  a  certain  common  value  of  k.  If  a  unique  curve  is 
obtained  in  this  way,  it  is  most  convenient  at  this  point  to 
try  to  reproduce  it  by  starting  from  the  theoretical  end, 
plotting  a.gainst  k  families  of  curves  of  /7"(K^)L^  (k)  for 
various  not  too  unreasonable  values  of  R,  and  for  the  various 
allowed  It  is  surprising  what  a  wide  variety  of  curves  can 
be  predicted  or  "postdicted”  in  this  way. 


Q-  Application  to  the  Present  Case 

The  experimental  angular  distributions  of  Ueil3on  et  al. 
(9)  for  the  reaction  33  11  (d,n)C/*  **(4.43  liev  level),  at 
bombarding  energies  of  -J-,  1,  1-J-,  and  2  Hev,  have  been  analysed 
in  this  way,  except  that  no  attempt  ha3  been  made  to  subtract 
3ome  homogeneous  background  from  each  curve.  (These  results 
agree  with  those  obtained  by  other  groups,  and  which  have  been 
considered  as  one  of  the  central  pieces  of  evidence  for  the 
heavy-particle  stripping  hypothesis.)  The  curves  obtained  as  a 
function  of  k  after  renormalization  are  shov/n  in  Rigure  6. 

It  turns  out  that  this  is  a  borderline  case  with  respect 
to  the  test  described  above.  Remembering  that  the  average 
ordinates  have  been  renormalized  ad  hoc,  so  that  only  the 
agreement  in  shape  and  in  abcissae  of  features  tests  the 
hypothesis,  their  coalescing  to  a  single  curve  is  not  too 
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convincing.  But  it  is  not  so  Bad  that  if  a  suitable  theoretical 
expression  Y/ere  put  forward  it  could  not  conceivably  be 
considered  as  agreeing  approximately  with  all  four  of  them. 

Inasmuch  as  they  oan  be  described  as  coalescing  to  a 
single  curve,  this  curve  is  described  as  follows:  a  smooth 
monotonic  decline  from  the  smallest  k  reached  (k  ~.27), 
levelling  off  to  a  flat  minimum  at  lc  -.80  to  .84;  then  a 
gradual  rise  from  k  —.84  to  k  ~  *91 .  This  is  what  must  be 
predicted  theoretically  if  one  is  to  do  without  other  processes 
than  deuteron  stripping. 

Curves  of  the  quantities  If  (K^ ) j/  V)  and  7T(K )  j  ?X( /^  )  for 
radii  H  of  2,  3,  4,  5,  6,  7,  8,  and  9  fermis  have  been  plotted 
as  functions  of  k,  to  decide  whether  an  unexpected  value  of  R 
or  of  the  ratio  P7/P,  might  produce  anything  like  the  average 
experimental  curve  described  above. 

PolloY/ing  the  proposal  made  in  Section  I)  for  introducing 
Coulomb  effects  into  (d,n)  stripping  theory,  the  quantities 
IT  (Kb)ff,‘(  f  )/{  1  and  )sf(  {)/{'  have  also  been  plotted 

for  the  same  wide  choice  of  radii.  Some  of  these  graphs  are 
presented  in  figures  7  to  10. 

Detailed  examination  of  all  these  curves  lias  revealed  that 
neither  the  pure  stripping  nor  the  Coulomb  function  curves,  for 
any  ratio  P^/EJ  ,  can  reproduce  the  general  features  of  the 
aggregate  experimental  curve  described  above.  Those  few  curves 
which  show  the  required  rise  from  k  -.30  to  .97  all  show  a 
marked  rise  from  k  -  .27  to  a  peak  at  some  higher  k,  in 
contrast  to  the  required  smooth  decline  from  k  =  .27  to  .80. 


m 

* 

* 


73 


•3 


.3 


74 


75 


•  3 


X  10  Cm 


4 


5 


.8 


76 


77 


Thus  it  seems  that  there  really  is  something  different 
from  deuteron  stripping,  with  or  without  Coulomb  corrections, 
talcing  place  here.  Besides  heavy -particle  stripping,  a 
combination  of  compound  nucleus  effects  with  deuteron  stripping 
remains  a  possibility. 


!•;-  Conclusions  about  Coulomb  Effects 

The  method  suggested  for  the  inclusion  of  Coulomb  effects 
in  (d,n)  stripping  has  therefore  not  explained  the  anomalous 
distributions  of  the  B  ^ (d,n)  reaction,  nevertheless  it  remains 
of  some  interest  in  itself. 

The  modifications  introduced  into  the  angular 
distributions  by  this  procedure  can  be  seen  in  the  accompanying 
graphs.  The  functions  */V  start  from  zero  at  (*-  0,  rise  to 
a  first  peak,  decrease  to  zero,  rise  again  and.  again  to  smaller 
and  smaller  peaks,  jU3t  a3  the  functions  V  1  do,  but  they 
reach  each  peak  later,  i.e.  at  a  higher  -value,  than  the 
corresponding  W  ) ;  and  because  of  the  factor  p  in  the 
denominator,  the  peaks  are  neither  as  high  nor  as  sharp;  the 
area  under  the  curve  is  pushed  out  to  greater  C •  This  peak 
shift  is  larger  for  larger  ? ,  the  Coulomb  functions  becoming 
more  and  more  different  from  their  °l  -  0  limiting  case  the 
spherical  Be33el  functions.  Remembering  that  at  a  given 
bombarding  energy  the  physical  range  of  0  from  0°  to  130* 
restricts  us  to  looking  at  a  fixed  corresponding  range  of  k 
(e.g.  from  k  *  .27  to  k  =  1.01  for  the  B  ,f  reaction  at  2  Rev; 
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see  Figure  6),  vie  can  readily  describe  qualitatively  the 
general  effects  of  these  modifications  on  angular  distributions: 
decrease  in  total  cross-section,  decrease  in  sharpness  of  the 
pattern*s  features,  shift  of  the  pattern  towards  higher  9  . 

This  shift  in  pattern  is  obviously  more  considerable  for 
larger  7,  hence  for  larger  Z.  It  can  also  be  seen  that  it 
decreases  with  increasing  bombarding  energy.  For  one  tiling  7 
is  inversely  proportional  to  k.  But  more  important  is  the  fact 
that  it  is  a  unique  shift ,  in  the  curve  of  ^  as  a  function  of 
k,  corresponding  to  a  fixed  k-interval;  but  a  fixed  k-interval 
corresponds  to  a  smaller  and  smaller  0 -interval  as  the 
bombarding  energy  increases,  as  a  consequence  of  Equation  PI. 
Hence  this  shift  in  the  pattern  behaves  quite  reasonably. 

Ih  fact,  these  modifications  are  just  what  is  expected  on 
the  basis  of  qualitative  arguments  already  presented,  and  what 
is  found  by  vastly  more  complicated  calculations  for  (d,p) 
reactions • 

It  should  therefore  be  interesting  to  test  the  predictions 
of  this  easy-to-apply  partial  correction  (considering  that  some 
tables  and  graphs  of  f )  are  now  available,  (32,  40)  ) 
against  the  more  elaborate  and  diff i cult -to- compute  Coulomb 
formulas  developed  for  (d,p)  reactions,  and  against  experiment 
in  cases  where,  unlike  the  case  presently  studied,  pure 
deuteron  stripping  still  3eems  to  prevail,  modified  by  Coulomb 
forces  but  not  spoiled  beyond  recognition  by  other  processes. 
Hot  only  because  of  the  enormous  advantage  in  computational 
simplicity  of  course,  but  because,  for  (d,n)  reactions,  the 
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methods  which  consider  Coulomb  effects 

/ 

whole  are  certainly  open  to  criticism, 
beginning  of  Section  IT. 


on  the  deuteron  as  a 

4 

as  pointed  out  in  the 


3-  The  Compound  nucleus  process 

There  remains  the  possibility  of  a  significant 
contribution  from  processes  involving  compound  nucleus 

U  /x 

formation,  in  the  B  (d,n)C  reaction.  This  has  always  been 
considered  unlikely  because  there  are  no  levels  of  C 13  in  the 
energy  region  reached  by  these  experiments.  (The  nearest  are 
two  levels  of  unknown  properties  at  excitations  corresponding 
to  2.18  and  3.08  Kev  bombarding  energy  for  the  B  "  +  d  reaction 
(Ajzehberg-Selove  and  Lauritsen,  41).  )  Bor  the  same  reason  it 
would  be  difficult  to  calculate  the  angular  distributions  and 
correlations  to  be  expected  on  this  basis,  since  we  are  neither 
on  a  resonance  nor  in  a  region  where  many  closely-spaced  levels 
contribute  statistically.  The  only  information  that  seems  well 
established  in  such  a  situation  is  that  angular  distributions 
should  be  symmetrical^ about  3  -  90°  (Wolfenstein,  42).  It  is 
al30  generally  accepted  that  the  distributions  should  not 
exhibit  features  as  sharp  as  those  characteristic  of  direct 
reactions • 

Y/ard  and  Grant  (43),  having  obtained  at  a  bombarding  energy 

of  .6  2Iev  an  angular  distribution  that  looked  roughly 

symmetrical  about  90°,  had  interpreted  it  by  compound  nucleus 
•f*  see  footnote,  page  81. 
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theory,  with  reservations,  ilore  accurate  and  extensive 
measurements  have  since  invalidated  this. 

A  general  treatment  should  include  both  compound  nucleus 
and  direct  interaction  contributions,  possibly  with 
interference  betv/een  them.  Thomas  (45)  ha3  begun  to  explore  the 
theory  of  such  messy  situations,  but  this  necessarily  involves 
even  more  loose  ends  than  the  work  reviewed  in  Chapter  2, 
especially  since  the  conditions  are  not  favorable  to  an 
unambiguous  treatment  of  either  of  these  two  processes. 

If  it  were  simply  a  matter  of  adding  without  interference 
a  deuteron  stripping  cross-section  and  a  cross-section 

p 

symmetrical  about  90  but  otherwise  undetermined,  one  might 
try  a  very  rough  analysis  of  the  experimental  angular 
distributions  in  the  following  way.  Starting  with  a  smoothed- 
out  version  of  the  experimental  curves,  attribute  the  whole 
&  90°  part  of  the  distribution  to  compound  nucleus  formation; 

a  corresponding  part  of  the  distribution  for  0  <C90°  can  then 
be  similarly  assigned  to  compound  nucleus  formation  because  of 

o 

the  rule  concerning  symmetry  about  90  ;  finally  the  remainder 

of  the  distribution  for  Q  <  90°  can  be  checked  for  consistency 

with  deuteron  stripping,  and  the  part  assigned  to  compound 

nucleus  formation  examined  to  determine  its  plausibility. 

Unfortunately  even  this  much  over-simplified  program  can 

not  be  carried  out  because  data  are  not  available  for  d  y  145°. 

Hence  the  program  fails  also  for  Q  <  35&,  and  it  is  not 

possible  to  decide  whether  a  segment  of  curve  ranging  only  from 
6  o 

35  to  90  is  part  of  a  stripping  curve  or  not.  lloreover  the 
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curve  for  llev  ‘bombarding  energy  shows  a  minimum  around  60 & 

0  O 

and  a  rise  from  60  to  100  ,  which  can  not  be  dealt  with  on 
these  assumptions.  Ho  conclusions  can  be  drawn  from  this 
concerning  the  importance  of  compound  nucleus  effects. 

Qualitatively,  an  explanation  along  those  lines  would 
require  such  effects  to  be  of  the  same  order  of  magnitude  a3 
the  deuteron  stripping  contribution  or  larger.  Ignoring 
interference}*  the  predicted  compound  nucleus  angular  distribution 
would  have  to  be  fairly  isotropic  for  bombarding  energies  of  -J- 
and  1  LZev;  and  to  show  a  minimum  at  90*  with  a  gradual  rise  on 
each  side  to  about  li  or  1-J-  times  the  mi ni mum  ordinate  at 
0^50  and  150°  ,  for  bombarding  energies  of  l-J-  and  2  Hev. 

1”  In  some  cases  there  may  be  terms  antisymmetric  about  9=-  90°, 
due  to  interference  between  two  levels  of  opposite  parity.  This 
possibility  is  not  included  in  the  rough  treatment  of  this 
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